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SUMMARY 


Previous  experimental  and  theoretical  work  on  atoms  in  solid  hydrogen  (Hj)  leads 
to  the  following  conclusions:  (1)  large  concentrations  of  atoms  can  exist  in  molecular 
solids  at  low  temperatures;  (2)  it  is  easier  to  cool  samples  with  known  radioactivity  than 
might  be  imagined;  (3)  the  atoms  are  inhomogeneously  distributed  throughout  the 
sample;  and  (4)  application  of  a  magnetic  field  may  increase  the  concentration  of  atoms 
that  can  be  stored  at  a  given  temperature  or  their  storage  time. 

Recent  work  has  focused  on  maximizing  the  stored  energy  in  hosts.  These 
hosts  were  divided  into  two  categories  depending  on  the  method  used  to  form  the  crystal: 
crystallization  ft:om  the  liquid  phase,  c-Hj,  or  deposition  directly  from  the  vapor  onto  a 
cold  substrate,  a-H2.  In  both  categories  of  solid  H2,  the  concentration  of  metastable 
particles  (atoms)  depends  on  the  diffusion  coefficient  since  it  controls  the  recombination 
rate.  Work  done  under  this  project  can  be  divided  into  four  areas:  (1)  quantification  of  the 
diffusion  and  recombination  coefficients  for  three-dimensional  (3-D)  hydrogen  hosts  of 
both  C-H2  and  a-Hj  types;  (2)  development  of  a  new  model  of  the  recombination  process 
that  predicts  observed  but  previously  unexplained  atom  growth  curves  and  suggests  a  new 
approach  for  obtaining  high  atom  densities;  (3)  development  of  a  new  model  for  the  in 
situ  production  of  atoms  in  a  solid  hydrogen  host;  and  (4)  new  nuclear  magnetic 
resonance  (NMR)  measurements  on  hydrogen  monolayers  on  activated  carbon  fibers  that 
yield  the  molecular  diffusion  coefficient  for  two-dimensional  (2-D)  hydrogen  and  an 
enhanced  ortho-para  conversion  rate. 

There  are  three  fertile  areas  for  future  investigation.  The  diffusion  coefficients 
obtained  in  3-D  hydrogen  hosts  were  for  temperatures  just  below  the  triple  point 
temperature  and  for  solids  crystallized  from  the  liquid  state.  Extension  of  these 
measurements  on  vapor  deposited  solids,  from  the  triple  point  to  the  lowest  possible 
temperatures,  would  be  very  important  to  the  High  Energy  Density  Materials  (HEDM) 
program.  In  addition,  use  of  non-conventional  NMR  techniques  to  extend  both  sets  of 
measurements  to  the  5  K  vicinity  would  determine  the  molecular  hopping  frequencies 
unambiguously.  The  values  of  the  diffusion  coefficient  obtained  for  2-D  hydrogen  were  so 
large  that  it  is  impractical  to  think  of  storing  metastables  in  surface  layers.  These  values 
should  be  checked  because  two  distinct  motional  processes  can  change  the  width  of  NMR 
lines;  (1)  actual  hopping  or  (2)  motion  of  the  spins  without  hopping.  In  the  latter  case, 
the  motional  narrowing  of  the  NMR  lines  does  not  give  information  regarding  hopping 
and  no  conclusion  regarding  the  storage  of  metastables  in  surface  layers  could  be  drawn 
from  the  data. 
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1.0  INTRODUCTION 


The  density  of  energy  that  can  be  stored  in  a  host  hydrogen  solid  depends  on  the 
diffusion  of  the  metastable  particles  that  carry  the  energy  and  that  of  the  host  solid.  Thus 
diffusion  of  metastables  in  H2  and  that  of  the  host  solid  molecules  is  of  critical  interest  to 
the  High  Energy  Density  Materials  (HEDM)  program. 

The  work  in  the  HEDM  program  described  in  this  report  concentrated  on  four 
aspects  of  high  energy  density  materials:  (1)  diffusion  of  molecules  and  atoms  in  3-D 
hydrogen  hosts;  (2)  suppression  of  atomic  recombination  in  H2  hosts;  (3)  in-situ 
production  of  atoms  in  H2  hosts;  and  (4)  diffusion  of  H2  molecules  in  2-D  monolayers. 

The  primary  technique  for  studying  H2  hosts  has  been  nuclear  magnetic  resonance 
(NMR),  one  of  the  few  experimental  methods  that  can  study  diffusion  microscopically 
since  the  hopping  motion  of  particles  with  nuclear  spin  causes  a  predictable  dephasing  of 
the  NMR  signal  so  that  the  diffusion  coefficient  can  be  measured  in  situ,  negating  the 
need  for  a  radioactive  tracer  element.  In  addition,  cracks  (or  other  extended  defects)  can 
actually  destroy  radioactive  tracer  experiments  but  have  no  effect  on  NMR  experiments. 
This  is  true  when  there  is  only  a  small  fraction  of  the  diffusing  atoms  (molecules)  in  the 
cracks  at  any  one  time  since  all  atoms  have  equal  weight  in  the  NMR  signal,  independent 

of  their  speed.  . 

The  solid  molecular  hydrogens  represent  an  especially  interesting  class  of  solids  in 

which  to  study  diffusion  since  there  are  three  isotopes  (H,D,  and  T)  from  which 
molecules  with  masses  of  2, 3,  4 , 5,  and  6  can  be  formed.  The  molecules  are  electrically 
identical  so  that  the  intermolecular  interactions  between  the  various  types  of  molecules 
are  the  same.  Since  the  molecular  masses  differ,  the  zero  point  effects  are  different,  the 
triple  point  of  H2  being  14  K  while  that  of  T2  is  20.6  K,  leading  to  very  different  effective 
potentials.  To  date,  there  have  been  several  experimental  investigations  pertinent  to  this 
work  [1-7],  but  only  one  theoretical  investigation  [8]  of  the  molecular  self-diffusion 
coefficient  in  the  solid  hydrogens.  New  data  on  molecular  diffusion  is  presented  and  new 
theoretical  models  are  used  to  extract  the  atom  diffusion  coefficient  so  that  it  can  be 
compared  with  the  molecular  diffusion  coefficient. 

1.1  Diffusion  in  3-D  Hydrogen  Hosts 

The  new  NMR  data  on  the  molecular  diffusion  coefficient  was  taken  on  the  solids 
HD,  D-T,  and  Tj.  This  data,  combined  with  existing  data  on  H2  [1-4],  HD  [1]  and  D2  [5- 
7],  yields  a  more  complete  picture  of  diffusion  in  the  solid  hydrogens.  The  complete  set  of 
data  is  used  empirically  to  obtain  the  activation  energy  as  a  function  of  isotope  mass.  A 
new  method  is  proposed  for  the  temperature  scaling  of  the  diffusion  data  from  one  isotope 
to  predict  the  behavior  of  another. 

Comparison  of  the  diffusion  coefficients  of  molecules  with  those  of  atoms  [9-11]  in 
sohd  hydrogen,  deduced  from  the  data  on  atomic  recombination  and  obtained  from  the 
analysis  here  and  earlier  by  Souers  [12]  shows  that  the  diffusion  coefficients  for  atoms  and 
molecules  are  comparable  in  the  solid  hydrogens.  This  is  to  be  expected  if  the  diffusion  of 
both  atoms  and  molecules  is  dominated  by  vacancy  exchange.  Finally,  the  scaling  used 
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for  the  molecular  diffusion  coefficient  is  used  for  the  atom  diffusion  coefficient  to  predict 
the  diffusion  coefficient  of  T  atoms  in  solid  Tj.  Calculation  of  the  atomic  recombination 
coefficient  in  Tj  enables  an  estimate  of  the  steady-state  atom  population  to  be  made.  This 
estimate  is  in  good  agreement  with  that  inferred  from  analysis  of  the  ortho-para 
conversion  data  [13]  in  solid  Tj. 

The  diffusion  coefficient  obtained  here  is  not  a  directly  measured  quantity  but  is 
deduced  from  measurements  of  the  damping  of  NMR  spin-echo  amplitudes  for  molecules 
and  the  electron  spin  resonance  (ESR)  signal  loss  in  time  (due  to  atomic  recombination) 
for  atoms.  We  will  review  the  expressions  that  will  be  used  to  obtain  the  diffusion 
coefficients  from  the  data  and  the  Ebner  and  Sung  theory  [8]  of  molecular  diffusion  in  Hj. 

1.2  Suppression  of  Atom  Recombination 

In  analyzing  the  present  results  on  diffusion,  it  was  found  that  both  atom  diffusion 
and  molecule  diffusion  in  the  solid  hydrogens  were  dominated  by  vacancies.  Thus  when 
an  atom  (or  molecule)  moves  by  hopping  into  an  available  vacancy,  it  creates  another 
vacancy  at  the  site  just  left.  Using  the  dependence  of  the  diffusion  coefficient,  D(T),  on 
the  concentration  of  vacancies,  it  was  concluded  that  as  atoms  were  created  from 
molecules  by  the  beta  decay,  the  available  vacancies  would  start  to  fill  up  since  the 
destruction  of  one  molecule  creates  one  vacancy  and  two  atoms.  A  quasi-equilibrium 
theory  was  obtained  to  predict  the  number  of  vacancies  present  in  solid  hydrogen  at  any 
temperature  as  atoms  are  created  from  molecules. 

1.3  Production  of  Atoms  In  Situ 

Hydrogen  solid  hosts  containing  tritium  such  as  Ta,  DT,  HT,  or  even  Hj,  HD,  and 
Dj  doped  with  some  amount  of  Tj  display  very  interesting  thermal  [14]  responses  such  as 
spontaneous  and  triggered  energy  releases  and  optical  properties  [15]  such  as  continuous 
emission  and  triggered  optical  pulses  at  low  temperatures.  It  is  known  that  these  solids 
can  store  hydrogen  atoms  but  how  many  of  these  metastable  excitations  can  exist  at  any 
given  temperature  is  still  unknown.  The  atoms  stored  have  their  origin  in  the  beta  decay 
of  the  tritium  nucleus.  In  the  gas  phase,  one  beta  decay  produces  about  800  atoms  [15] 
but  this  number  could  be  different  in  the  solid. 

Experiments  on  such  solids  have  revealed  many  interesting  phenomena.  Also,  the 
solids  themselves  are  intrinsically  interesting  because  they  represent  potential  energy 
storage  systems,  possibly  for  space  propulsion,  and  also  are  the  designated  cryogenic  fuel 
for  inertial  confinement  fusion  experiments.  In  spite  of  the  great  amount  of  effort  spent 
on  some  of  the  quantum  aspects  of  Hj,  the  materials  properties  of  Hj  are  not  particularly 
well  known  when  compared  to  common  substances.  To  date,  there  has  been  little 
theoretical  development  of  the  subject  of  this  report,  namely  the  production  of  atoms  in 
solid  hydrogen  due  to  tritium  beta  decay  although  Rosen  [16]  and  Zeleznik  [17]  did  model 
the  unusual  thermal  responses  observed  by  Webeler  [18]  on  Hj  containing  1%  (or  less)  Tj. 
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1.4  Diffusion  in  2-D  Monola)/ers 

Pulsed  NMR  techniques  were  used  to  study  H2  and  HD  molecular  monolayers 
absorbed  on  activated  carbon  fibers  (ACF)  over  the  temperature  range  from  4.6  K  to  20  K. 
These  fibers  have  an  effective  absorption  area  of  2000  m^  per  gram  so  that  for  monolayer 
coverage,  10'®  molecules  per  m^,  1  gram  of  ACF  holds  about  1  liter  at  standard 
temperature  and  pressure  (STP)  of  hydrogen  gas.  The  Fourier  transform  of  the  NMR 
lineshapes,  the  free  induction  decay  (FID),  and  the  measured  height  of  the  FID 
immediately  after  the  90°  pulse  (a  quantity  that  is  proportional  to  the  number  of  spins  in 
the  sample)  were  recorded  as  functions  of  time  and  temperature.  The  time  dependence  of 
the  NMR  signal  in  makes  it  possible  to  measure  the  ortho  (spin  1)  to  para  (spin  0) 
conversion  time  since  only  the  ortho  molecules  have  an  NMR  signal.  The  FID  gives  the 
lineshape  and  enables  one  to  study  the  spin  dynamics  of  the  molecules.  Analysis  of  the 
NMR  data  is  based  on  comparisons  with  that  of  0H2-PH2  mixtures  of  bulk  3-D  solid 
hydrogen. 
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2.0  TECHNICAL  DEVELOPMENTS 


Technical  Developments  is  divided  into  four  sections:  (1)  diffusion  measurements 
in  3-D  hydrogen  hosts;  (2)  suppression  of  recombination  of  atoms  in  3-D  hydrogen  hosts; 
(3)  production  of  atoms  in-situ;  and  (4)  monolayers  absorbed  on  activated  carbon 
fibers.  In  all  sections,  theoretical  advances  are  included  along  with  the  experimental 
work. 

2.1.  Diffusion  Measurements  in  3-D  Hydrogen  Hosts 

To  correctly  interpret  the  NMR  data,  it  was  necessary  to  develop  a  new  consistent 
model  of  the  hop  motion  of  atom  and  molecules  in  hydrogen  hosts.  This  model  is 
described  below  and  then  used  to  analyze  the  data  on  3-D  diffusion  and  later  the  data  on 
diffusion  in  monolayers  of  H2  absorbed  on  activated  carbon  fibers. 

2.1.1.  Theoretical  Background.  Consider  only  the  monovacancy  limit  where  the 
concentration  of  vacancies,  c^,  is  very  small.  Let  To  be  the  basic  rate.  It  is  the  probability 
per  unit  time  that  a  host  molecule  (e.g.,  Hj)  will  hop  to  a  vacant  neighboring  site.  Then 
the  vacancy  diffusion  coefficient,  D^,  the  diffusion  constant  associated  with  vacancy 
motion,  is: 

(1) 

where  Z  is  the  number  of  neighbors,  a  is  the  nearest  neighbor  distance,  taken  to  be  the 
distance  of  one  hop,  and  3  is  the  number  of  dimensions.  For  hexagonal  close  packed  (hep) 
and  face  centered  cubic  (fee)  lattices,  Z  =  12.  The  probability  per  unit  time  that  a 
vacancy  will  hop  is  ZTo  because  a  vacancy  has  Z  available  sites.  On  average,  one-sixth  of 
the  jumps  take  the  vacancy  in  a  given  direction. 

This  will  be  applied  to  the  diffusion  of  molecules.  The  probability  that  any  site 
neighboring  a  given  molecule  is  vacant  (and  a  molecule  can  hop  to  it)  is  c^.  Thus 
approximately,  the  diffusion  coefficient  for  a  molecule  is  cT^^,  though  in  fact  this  is  not 
exactly  correct  since  the  molecule  under  consideration  and  the  vacancy  can  undergo 
multiple  scattering.  (For  example,  the  molecule  and  vacancy  could  exchange  places  twice, 
which  would  lead  to  no  effect.)  The  diffusion  coefficient  for  a  molecule  is  written  as 

(2)  D, 

where  f„  is  the  (tracer)  correlation  factor.  It  can  be  shown  that 


and  f,„  =  0.781  for  an  fee  lattice  and  certain  motions  on  an  hep  lattice. 
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The  molecular  hop  rate,  r„,  is  defined  as  r„  =  c^ZT^  =  1/x^  where  is  the 
correlation  time.  Except  for  details  like  one  can  think  of  as  the  effective  hop  rate  of 
a  host  molecule  to  any  site  and 


(4) 


= 


/  '7\ 


(6 


T 


It  is  the  hop  rate,  r„j,  that  Ebner  and  Sung  calculated  [8]. 

Since  the  data  do  not  yield  D  directly,  it  is  important  to  relate  D  to  the  measured 
quantities.  The  NMR  data  consist  of  measurements  of  the  Fourier  transform  of  the  NMR 
lineshapes,  the  free  induction  decay,  the  measured  height  of  the  FID  immediately  after  the 
90°  pulse  (a  quantity  that  is  proportional  to  the  number  of  spins  in  the  sample),  and  the 
temperature  dependence  of  the  NMR  time  constant  (the  transverse  relaxation  time) 
denoted  Tjj.  This  quantity  is  usually  denoted  by  the  s3mibol  Tj  but  in  this  report  Tjj  is 
used  to  avoid  confusion  with  molecular  tritium  in  this  section.  It  is  roughly  equal  to  the 
inverse  linewidth,  from  spin-echo  experiments.  The  measured  quantity  T^^  is  related  to 
the  correlation  time,  by  the  expression  [19]  for  fee  and  hep  lattices: 


(5) 


M, 


8.67 


hx. 


where  M2  is  the  rigid  lattice  second  moment  and  the  pure  number,  b,  depends  on  the 
product  of  the  Larmor  frequency  of  the  experiment,  S),  and  the  correlation  time  x^with: 
b  =  38.3,  when  S)x,  <  <  1  and  b  =  11.6,  when  B)x,  >  >  1.  All  the  data  taken  here  are  in 
the  latter  region,  so  b  =  11.6.  Thus  =  1.338  M2X(..  Including  the  factor  f„,  the 
diffusion  coefficient  for  molecules,  D,„,  will  be  calculated  from  the  expression: 


(6) 


Z)„  =  [1.338 1  yJ  = 


The  net  effect  of  the  two  corrections  on  D„,  thus  amounts  to  less  than  5%  when  compared 
to  expressions  used  previously  in  the  literature.  The  actual  lattice  spacings  for  the 
isotopes  [12]  are  used  to  calculate  from  the  data. 

The  single  theoretical  study  pertinent  to  this  work  was  performed  by  Ebner  and 
Sung  [8].  The  Ebner  and  Sung  paper  was  reviewed  and  augmented  in  the  paper  by  Rail,  et 
al.  [4].  Ebner  and  Sung  considered  two  possible  mechanisms:  (1)  thermally  activated 
diffusion  over  a  barrier  into  a  vacancy  and  (2)  tunneling  through  the  barrier  into  a 
vacancy.  The  former  process  would  dominate  at  high  temperatures  while  the  latter 
process  would  dominate  at  low  temperatures.  A  diffusion  coefficient  of  the  form 


(7) 


D  =  Dq  exp  (-  •^) 
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was  predicted  in  both  temperature  regions,  with  the  numerical  values  of  both  Dq  and 
depending  on  the  specific  Emit.  Although  Ebner  and  Sung  did  not  consider  the  hop  of 
atoms  in  solid  Hj,  their  results  for  molecules  in  were: 


£)  =  (2  X  10-^  c/w^/s)  expf- 

\  T  > 


[/ow  7] 


where  m  is  the  mass,  E^  is  the  energy  to  create  a  vacancy,  E,,  is  the  barrier  height,  E^  is  the 
energy  of  the  first  tunneling  level,  and  Jj  is  the  splitting  of  the  first  tunneling  level. 


2.\2  The  Diffusion  of  Atoms.  There  have  been  no  direct  measurements  of  D  for 
atoms.  However  measurements  of  the  recombination  coefficient  can  be  analyzed  to  yield 
the  atom  hop  rate  and,  hence,  the  diffusion  coefficient.  Defining  2  a  as  the  atom 
recombination  coefficient,  Vq,  the  volume  per  host  molecule  (aVV’2),  Z  the  number  of 
nearest  neighbors,  the  time  to  hop  to  a  specific  site,  t,  the  recombination  time,  and  a 

constant  of  order  unity,  the  recombination  coefficient  can  be  shown  to  be  equal  to: 


(10) 


2a 


If  it  is  assumed  that  >  >  Zt„  then  2a  =  2voZra„y,  where  =  Z/thy  is  the 
rate  for  an  atom  to  hop  to  any  site  and  thus  Z  times  larger  than  the  rate  to  hop  to  a 
specific  site.  The  factor  7?  is  obtained  because  there  are  Z  nearest  neighbors  to  an  atom 
and  Z  sites  leading  to  each  of  them.  Thus  from  the  recombination  coefficient,  the  atom 
hop  rate  is  obtained,  which  can  then  be  compared  to  the  molecular  hop  rate. 

2.1.3  Comparison  with  Experiment.  At  the  time  the  Ebner  and  Sung  calculation 
was  completed,  NMR  data  on  the  diffusion  coefficient  of  Hj,  HD,  and  Dj  existed. 
Comparison  of  the  theoretical  expressions  for  Hj  and  Dj  with  existing  experimental  data 
gave  very  good  agreement  in  the  high  temperature  regime  but  no  crossover  into  the  low 
temperature  regime  was  observed  in  the  data.  Ebner  and  Sung  based  their  comparison 
with  experiment  on  an  analysis  of  Bloom's  NMR  data  by  Moriya  and  Motizuki  [20]. 

Although  the  details  of  their  analysis  differ  from  those  in  this  report,  Ebner  and 
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Sung  obtained  the  diffusion  coefficient  for  Hj,  D  =  (2  x  10'®)exp(-190/T)  cmVs,  a  result 
that  agreed  very  well  with  their  prediction.  The  activation  energy  of  288  K  had  been 
obtained  for  Dj  from  the  work  of  H.  Meyer  and  co-workers  [5]  and  agreed  well  with  that 
calculated  by  Ebner  and  Sung.  Thus  the  activation  energies  for  both  Hj  and  Dj,  obtained 
from  experiment,  agreed  extremely  well  with  the  calculations. 

Ebner  and  Sung  did  not  compare  their  calculation  with  the  existing  data  on  HD 
[1] ;  possibly  no  calculation  was  done  for  HD.  Considering  how  different  the  existing  HD 
data  was  from  that  of  Hj  and  Dj,  it  is  unlikely  that  there  could  have  been  any  agreement. 
Both  the  activation  energy  (EJ  and  the  prefactor  Dq  in  HD  appear  to  be  anomalous  when 
compared  to  the  Hj  and  Dj  values.  Bloom  comments  that  there  were  difficulties  in 
establishing  good  thermal  contact  with  the  HD  sample  that  did  not  exist  for  Hj,  probably 
leading  to  an  erroneous  determination  of  the  activation  energy  (and  hence  Dq)  for  HD. 
Since  the  new  experimental  values  for  Ej  and  Dgin  HD  do  not  appear  to  be  anomalous 
when  compared  to  Hj  and  Dj,  they  were  used  and  the  earlier  values  were  ignored. 

The  agreement  of  the  Ebner  and  Sung  theory  with  activation  energies  and  even 
prefactors  obtained  from  the  existing  data  on  Hj  and  Dj  appeared  excellent  since  the 
calculation  involved  no  adjustable  parameters.  However,  if  the  theoretical  result  is  used 
to  generate  numerical  values  of  D  in  the  temperature  interval  from  14  K  to  1 1  K,  the  range 
of  the  H2  experiments,  and  these  points  are  then  fitted  to  Equation  1,  one  finds  Dq  =  2.8  x 
10'®  cmVs  and  E^  =  139  K,  values  that  do  not  agree  with  experiment.  This  points  to  a 
particular  problem  with  the  calculation,  namely  that  a  crossover  to  quantum  tunneling  is 
predicted  about  14.7  K,  a  temperature  just  above  the  triple  point  of  Hg-  This  crossover  has 
never  been  observed  experimentally  in  any  work  on  molecules.  Since  this  prediction  has 
never  been  verified,  the  quantum  part  of  the  Ebner  and  Sung  calculation  may  be 
incomplete. 

2.1.4  Experimental  Data.  New  measurements  of  the  molecular  self  diffusion 
coefficient  in  the  solid  hydrogens  have  been  made  by  conventional  pulsed  NMR 
techniques  in  Tj,  DT,  and  HD  and  the  original  measurements  in  Hj  have  been  verified. 

The  new  measurements,  at  30  MHz,  when  combined  with  earlier  ones  provide  a  more 
complete  picture  of  molecular  diffusion  in  these  quantum  solids  and  can  be  interpreted  to 
give  some  additional  insight  into  the  relationship  between  atomic  diffusion  and  molecular 
diffusion.  The  sample  configuration,  NMR  spectrometer,  and  details  of  the  experimental 
procedure  have  been  previously  published  [13] . 

Experimentally,  the  temperature  dependence  of  the  transverse  relaxation  time,  T^^, 
in  the  region  where  motional  effects  dominate  the  observed  linewidth  was  measured.  The 
samples  are  loaded  into  the  cell  as  liquids.  The  temperature  is  then  reduced  in  increments 
of  about  0.5  K.  A  Carr-Purcell-Meiboom-Gill  pulse  sequence  was  used  to  measure  T^^  at 
the  higher  temperatures  but  a  single  90°  pulse  was  used  when  T^^  became  too  short  for  the 
longer  sequence  to  be  practical.  Quadrature  detection  was  used  so  that  the  signals  had  to 
be  phase  corrected  and  then  extrapolated  to  t  =  0. 

All  other  workers  who  measured  the  diffusion  coefficient  in  molecular  hydrogen 
have  used  the  above  approach,  without  quadrature  detection,  with  the  exception  of  the 
work  on  Dj  by  Meyer  and  co-workers  [6].  They  used  a  more  sophisticated  approach 
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involving  NMR  rotating  frame  measurements,  a  technique  that  is  well  grounded, 
theoretically.  By  comparing  their  own  conventional  measurements  with  the  rotating 
frame  measurements,  they  demonstrated  the  values  obtained  by  the  two  techniques 
agreed  (to  within  experimental  error),  justifying  the  usage  of  the  less  sophisticated 
experimental  approach  used  here. 

The  new  NMR  measurements  of  in  HD,  D-T  and  Tj  are  given  in  Figure  1. 
Points  taken  in  the  liquid  phase  are  included  in  this  figure  but  are  not  included  in  the 
analysis  of  the  solid  data.  The  raw  data  have  been  converted  to  calculated  diffusion 
coefficients,  D(T),  in  Figure  2,  making  it  possible  to  determine  the  activation  energies  for 
all  the  isotopes.  Data  taken  on  the  isotopes  and  Dj  by  H.  Meyer  and  co-workers  [5,6] 
are  included  for  comparison  and  completeness. 


Figure  1 

Tjg  Versus  the  Reciprocal  Sensor  Temperature,  T'* 

The  plot  s3mibols  are:  for  HD  [A];  D-T  [■□];  and  Tj  [Mnj. 
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Points  in  the  liquid  phase  are  shown  but  not  included  in  the  analysis.  A  temperature 
correction  must  be  applied  to  the  Tj  and  D-T  points  due  to  their  self-heating. 

10-8 


10-9 

'vi 

^  I 

E  10-1 

Q 

lo-i 


lo-i 


lO-'^ 

0.04  0.05  0.06  0.07  0.08  0.09  0.10 

1/T(K) 

Figure  2 

D  versus  the  Reciprocal  Sensor  Temperature,  T  ‘ 

The  plot  symbols  are:  Tj  [Mo],  D-T  [■],  Dj  [V],  HD  [A],  and  [vj. 

One  might  expect  that  since  the  parameters  that  characterize  the  6-12  potential  in 
all  the  isotopes  are  the  same,  the  diffusion  coefficients  should  have  universal  behavior,  but 
they  don't.  The  behavior,  while  not  universal,  is  smooth,  showing  only  simple  thermal 
activation.  Accordingly,  Equation  7  was  used  to  fit  the  data  and  extract  values  of  and 
Do. 

Most  workers  in  this  field  quote  Dq  with  an  error  of  ±  100%  because  of  the 
approximate  nature  of  the  expressions  used  to  calculate  Dq.  Moreover,  a  small  error  in  E, 
can  make  a  big  change  in  Dq.  There  are  additional  problems  in  radioactive  samples. 

There,  temperature  gradients  exist  in  the  solid  due  to  the  heating  resulting  from  the  beta 
decay  of  the  triton.  The  results  presented  here  are  corrected  for  these  gradients  and 
average  temperatures  are  quoted.  Still,  these  gradients  can  lead  to  an  error  of  three  to  five 
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in  Dq.  The  temperature  correction  is  discussed  in  detail. 

For  both  Tj  and  D-T,  the  sensor  temperature  is  not  necessarily  equal  to  the  average 
sample  temperature.  In  the  experiments,  the  volume  of  the  samples  W'ere:  0.157  cm^  for  Tj 
and  0.172  cm®  for  D-T.  These  numbers  are  based  on  use  of  200  cm®  of  gas  at  STP.  The 
sapphire  cell  (a  cylinder)  had  an  inside  diameter  of  0.762  cm,  making  the  thickness  0.366 
cm  for  Tj  and  0.377  cm  for  D-T.  Thus  the  sample  radius  is  approximately  equal  to  the 
thickness.  Heat  could  be  removed  from  the  bottom  or  radially,  or  both.  Unfortunately, 
hydrogen  samples  have  a  tendency  to  pull  away  from  walls  as  they  are  cooled  below  the 
triple  point  so  it  is  not  certain  that  the  sample  remained  in  contact  with  the  walls.  The 
temperature  correction  will  be  estimated  in  the  best  case,  where  the  sample  remains  in 
contact  with  the  walls  and  in  the  worst  case  where  the  sample  pulls  away  from  the  walls. 

Based  on  molar  volumes  of  18.75  x  10**^  m®/mol  for  Tj  and  19.29  x  10'*^  m®/mol  for 
D-T,  the  samples  generated  1.745  x  10'®  W  (Tj)  and  8.723  x  10'®  W  (D-T).  Using  the 
thermal  conductivity  values  from  Collins  et  al.  [17],  with  0.35  W/mK  for  Tg  and  0.30 
W/mK  for  D-T  as  representative  values,  the  average  temperatures  for  the  best  case  (both 
radial  and  vertical  heat  flow)  are  given  by: 

Tj:  Tjvg  =  Tj^sor  +  0.488  K 

D-T:  T,^  =  T_,-^  0.237  K 

For  the  worst  case,  T^^  =  +  HL®/3k,  where  H  represents  the  self-heating  and  k  is 

the  thermal  conductivity  (with  L  the  length)  giving 

^2=  Ta^  =  Tjansor  +  0-^^^  K 

D-T:  Ta^  =  T_,  +  0.343  K 

Using  these  corrections,  it  is  now  possible  to  obtain  the  activation  energies  and  the 
appropriate  error  bars.  Measurements  on  DT  between  the  temperatures,  21.6  K  and  16.1 
K,  5deld  Fa  =  372  K  and  Dq  =  0.34  cm®/s.  Measurements  on  Tj  between  the  temperatures, 
21.7  K  and  17.7  K,  yield  Ea  =  432  K  and  Dq  =  1.8  cm®/s.  The  activation  energies  for  Tj 
and  D-T  are  based  on  the  best  case  heat  removal.  If  error  bars  are  derived  by  using  the 
worst  case  scenario,  one  has  for  Tj,  Ea  =  432  K  ±20  K  and  D-T,  Ea  =  372  K  ±10  K. 

Additional  errors  come  from  the  measurement  error  in  the  quantity  T^^  used  in 
Equation  6  to  calculate  D.  Because  of  the  very  high  signal-to-noise  ratio  in  the 
experiments,  the  error  in  determining  T^^  is  less  than  ±5%.  Calculated  values  of  Mj  were 
used  in  Equation  6,  and  since  these  values  depend  on  the  nearest  neighbor  spacing  (a),  the 
value  of  D  depends  on  (a)"^.  The  error  in  this  quantity  was  neglected  since  the  molar 
volumes  are  known  to  be  better  than  1%.  In  addition,  the  errors  in  the  measurement  of 
temperature  where  it  is  measured  and  regulated  to  better  than  50  mK  at  20  K  have  been 
neglected.  Thus  the  values  of  D  calculated  from  Equation  6,  due  to  random  errors  in  T^^ 
are  uncertain  to  ±5%.  For  extreme  cases,  this  also  translates  into  an  error  in  the 
activation  energy  of  ±5%,  also.  Combining  the  two  sources  of  error  for  Tj  and  D-T  one 
obtains  for  Tj,  E^  =  432  K  ±42  K  and  D-T,  E,  =  372  K  ±29  K.  A  random  error  in  Ej  of 
±5%  can  introduce  an  error  in  Dq  of  a  factor  of  2  or  3.  Thus  the  values  of  Dq  obtained 
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here  are  only  order  of  magnitude  estimates.  Th^  are  listed  below  to  two  significant 
figures  to  make  it  possible  to  calculate  individual  values  of  D  to  within  10%. 

Measurements  on  H2  and  HD  were  repeated.  The  results  on  Hj  agreed  well  with 
the  earlier  data  of  Bloom  [1]  and  Weinhaus  and  Meyer  [3]  and  the  more  recent 
measurements  by  Rail  etal.  [4].  The  new  results  on  HD  disagree  completely  with  the 
early  results  of  Bloom  [1] .  For  measurements  between  T  =  12.5  K  and  T  =17.1  K,  the 
new  data  )delds  =  250  ±  12  K  and  Dq  =  5.4  x  10'®  cmVs  (correct  to  within  a  factor  of  2 
or  3).  The  results  obtained  here  and  selected  results  of  other  workers  are  given  in  Table  I 
[18]. 


Table  1.  The  Activation  Energies  and  Diffusion  Coefficient  Prefactors  for  Molecules 


Isotope 

Ea(K) 

H(K) 

Do(cmVs) 

A 

H3 

HDinHa 

197  [3] 

197  [4] 

92.6 

2.8  X  10-® 

5.7  X  10-" 

1.731 

HD 

250 

116.7 

5.4  X  10-" 

1.471 

D2 

300  [7] 

138.3 

1.0  X  10-2 

1.224 

D-T 

372 

150.3 

0.34 

1.111 

L _ 

432 

164.8 

1.8 

1.00 

For  reference,  the  heat  of  sublimation,  H(K),  and  the  quantum  parameter  (A),  defined  in 
Equation  11,  are  included. 

2.1.5  Analysis  and  Interpretation.  In  this  section  data  taken  on  molecules  and 
atoms  in  the  solid  hydrogens  will  be  analyzed  and  interpreted.  Some  of  the  samples  were 
made  by  slow  crystallization  from  the  liquid  phase  while  others  were  formed  by 
deposition  of  vapor  on  a  cold  substrate. 

In  attempting  to  compare  data  on  molecules  in  different  quantum  solids,  it  is 
customary  to  use  reduced  variables,  an  approach  similar  to  that  used  in  applying  the  law 
of  corresponding  states  to  the  behavior  of  different  gases  in  an  attempt  to  see  the 
universal  features  of  all  gases.  The  reduced  variables  in  this  case  are  formed  from 
dimensionless  combinations  obtained  from  the  two  Lennard-Jones  parameters,  e  and  a 
[15] .  Unfortunately,  these  parameters  are  approximately  the  same  for  all  the  hydrogen 
isotopes  so  that  one  prediction  of  this  approach  would  be  that  the  product  of  the  diffusion 
coefficient  (D)  and  the  square  root  of  the  mass  (V m)  would  be  the  same  for  all  the 
hydrogen  isotopes.  This  is  clearly  incorrect  for  the  data  in  Figure  2.  For  instance,  by 
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drawing  a  horizontal  line,  it  can  be  seen  that  the  same  values  of  D  occur  for  widely 
different  values  of  1/T  or,  for  a  given  value  of  T,  the  diffusion  coefficient  of  the  lightest 
isotope,  H2,  is  far  larger  than  that  of  the  heaviest  isotope,  Tj.  This  is  a  clear  quantum 
effect.  A  different  method  of  comparing  the  isotope  data  is  required. 

What  was  needed  was  a  non-trivial  temperature  scaling  so  that  it  could  be 
determined  how  the  parameters  Dq  and  vary  from  isotope  to  isotope,  even  though  the 
potentials  are  essentially  the  same,  except  for  the  equilibrium  spacings  of  the  molecules. 
One  method  that  works  is  temperature  scaling  based  on  the  quantum  parameter  [15]. 
The  basic  idea  is  shown  in  Figure  3  where  the  triple  point  temperature  (Tjp)  of  each 
isotope,  multiplied  by  the  quantum  parameter  for  that  isotope,  is  plotted  as  a  function  of 
the  quantum  parameter.  A,  where 


(11) 


A  = 


^  i  K-nyr  ^0* 

6  AW.  e)'” 


where  Nq  is  Avogadro^s  number,  h  is  Planck's  constant,  a  and  e  are  the  Lennard-Jones 
parameters,  and  and  M,,  are  molecular  weights  of  the  atoms  in  the  molecule  (for 
heteronuclear  molecules),  with  M  being  the  molecular  weight  of  the  molecule.  The 
product  ATfp  is  roughly  constant  with  a  mean  value  of  22.8  and  an  rms  deviation  of  6% 
even  though  the  triple  point  temperatures  range  from  13.8  for  to  20.6  for  Tj.  Since  the 
triple  point  temperature  (essentially  an  energy)  scales  as  1/A,  the  same  A  dependence  for 
the  activation  energy  is  expected.  The  product  E^  A/16  is  plotted  in  Figure  3.  The  product 
Ej  A  is  constant  to  within  10%. 
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A  (quantum  par.) 

Figure  3 

Ttp  A  and  E,  A  versus  A 


This  approach  to  scaling  is  illustrated  in  Figure  4,  where  the  diffusion  coefficients 
for  all  the  isotopes  are  plotted  as  a  function  of  [TA]  '  =  l/Tj,  where  is  the  scaled 
temperature.  The  data  for  all  the  isotopes  fall  on  essentially  parallel  lines  indicating  the 
temperature  scaling  is  adequate  but  the  prefactors,  Dq,  are  all  different  and  no  suitable 
method  for  scaling  them  has  been  found  although  the  non-radioactive  isotopes  appear  to 
be  regular  and  the  radioactive  isotopes  are  anomalous. 
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Figure  4 

D  versus  =  TA 

The  plot  symbols  are:  Tj  [Mn],  D-T  [■□],  Dj  [V],  HD  [A],  and  Hj  [v]. 


The  measured  activation  energy,  at  high  temperatures,  is  composed  of  three  terms, 
according  to  Ebner  and  Sung  [8],  -  Ej,  where  E^  is  the  vacancy  energy,  E^  is 

the  barrier  energy,  and  Ej  is  the  energy  of  the  first  tunneling  level.  By  suitable  modeling,  it 
is  possible  to  estimate  the  three  contributions  individually.  Ebner  and  Sung  calculated  the 
lowest  energy  level  in  the  potential  well  for  and  obtained  the  value  of  E,  =  34.1  K,  a 
value  they  also  used  for  Dj.  If  the  shape  of  the  well  doesn't  change  appreciably,  then  the 
energy  levels  for  the  other  isotopes  should  be  simply  obtained  from  the  hydrogen  value, 
namely  E^ «  A,  where  A  is  the  quantum  parameter  (since  the  harmonic  oscillator  energy 
depends  on  l/V m.  Thus  the  tunneling  energies  are:  Hj  (34.1  K);  HD  (29.0  K);  Dj  (24.1  K); 
D-T  (21.9  K);  and  Tj  (19.7  K). 

The  energy  to  form  a  vacancy  in  a  perfect,  non-quantum  solid  hydrogen  is  210  K 
[21]  and  quantum  effects  should  reduce  this  value  proportional  to  l/Vm;  hence,  a  power 
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series  in  the  quantum  parameter  is  suggested.  Keeping  only  the  first  three  terms, 

(12)  £,  =  210  K-  K,{A)  +  K^(Af 

Ebner  and  Sung  gave  the  values  of  112  K  and  132  K  for  and  D2  so  both  and  Kj  can  be 
determined.  This  jdelds  the  values,  E^  =  112  K,  121  K,  132  K,  137  K,  and  143  K  for  H2, 
EID,  D2,  D-T,  and  T2,  respectively.  There  is  no  information  about  the  values  for  HD,  D-T, 
or  T2  in  the  literature. 

For  Eb,  there  appears  to  be  no  theoretical  guideline  similar  to  the  one  used  for  E^. 
The  experimental  value  of  E^,  the  calculation  of  E^  and  the  fit  to  E^,  are  used  to  obtain 
values  of  Eb  =  Ej  +  Et  -  E^  for  all  the  isotopes.  Ebner  and  Sung  calculated  Eb  =  119  K  and 
192  K  for  H2  and  Dj.  The  values  obtained  in  this  manner  are:  Eb  =  119  K,  158  K,  192  K, 
257  K,  and  309  K  for  H2,  HD,  D2,  D-T,  and  T2,  respectively. 

From  the  determination  of  the  hop  frequencies  for  the  isotopes,  it  is  possible  to 
calculate  the  hop  frequency  of  a  vacancy  at  the  triple  point  for  comparison.  All  the  solid 
hydrogens  have  about  the  same  Debye  temperature,  approximately  100  K.  This  translates 
into  a  Debye  frequency  of  (QJ  =  1.31  x  10^^  s\  four  orders  of  magnitude  larger  than  the 
vacancy  hop  frequency  at  the  triple  point  temperature,  as  seen  below. 

To  extract  the  hop  frequency  of  vacancies,  at  the  triple  point,  for  H2  one  uses  the 
expression: 


(13)  =  12  cJo  =  12  r„exp 

Solving  for  Tq  using  values  of  E^  taken  from  best  fit  gives,  yields  at  the  triple  points: 

ro=  2.75  X  10^  s'^  for  H2  using  triple  point  temperature  of  13.8  K. 

Tq  =  1.08  X  10®  s'‘  for  HD  using  triple  point  temperature  of  16.6  K 
ro=  0.62  X  10®  s'‘  for  Dj  using  triple  point  temperature  of  18.7  K 
ro=  1.21  X  10®  s'^for  D-T  using  triple  point  temperature  of  19.8  K 
ro=  0.86  X  10®  s'‘  for  Tj  using  triple  point  temperature  of  20.6  K. 

Thus  to  within  a  factor  of  4,  the  vacancy  hop  rate,  at  the  triple  point  temperature,  is  the 
same  for  all  the  hydrogen  isotopes. 

The  hopping  of  atoms  controls  the  steady-state  atom  population  trapped  in  solid 
hydrogen.  This  has  never  been  studied  directly  through  a  measurement  of  the  atom 
diffusion  coefficient  but  there  are  two  other  approaches  that  have  probed  the  steady-state 
atom  population  and  the  atom  dynamics.  These  are  ESR  studies  and  measurements  (and 
modeling)  of  the  ortho  to  para  conversion  of  molecules  of  T2  in  solid  Tj.  In  this  case,  there 
have  been  measurements  on  crystalline  solids  and  also  measurements  on  amorphous 
solids. 

Leach  and  Fitzsimmons  [9]  produced  H  atoms  in  crystalline  H2  by  electron 
bombardment  and  studied  atom  recombination  after  the  electron  beam  was  turned  off. 
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Defining  n  to  be  the  atom  density,  after  the  beam  is  off, 

(14)  ^  =  -  {2ay 

dt  ^  ^ 

where  2a  is  the  atom  recombination  coefficient.  They  could  measure  it  directly  from  the 
time  decay  of  the  ESR  signal,  which  was  proportional  to  the  density  of  atoms,  n.  They 
found  that 


(15)  a  =  ttjjcxp 

where  the  activation  energy  Ej,  =  195  ±10  K,  remarkably  like  that  found  for  molecules 
in  solid  Hj.  The  relationship  between  the  recombination  coefficient  and  the  atom  hop 
rate  (F),  Equation  10,  can  be  used  to  extract  F  from  the  measured  values  of  a.  After 
setting  the  activation  energies  to  a  common  value,  it  is  found  that 

(F)”  =  (5.9  X  10^5  sec*>xp(-197/r) 

This  can  be  compared  with  the  molecular  hop  rate  obtained  from  the  NMR  data 

(F)®  =  (1.56  X  10‘"  sec-‘)exp(-197  K/T) 

showing  that  the  atom  and  molecule  hop  rates  are  almost  identical  in  Hj,  differing 
only  by  a  factor  of  4. 

Other  important  information  about  crystalline  solids  is  provided  by  the  data  taken 
using  NMR  techniques  on  the  ortho  to  para  conversion  in  solid  Tj  [13] .  There  have  been 
no  direct  measurements  of  the  atom  recombination  coefficient  in  solid  tritium.  Because 
the  ortho-para  conversion  process  is  catalyzed  by  the  presence  of  atoms  [13],  the  rate  is  a 
sensitive  function  of  the  steady-state  atom  population.  Using  this  fact,  one  can  extract 
the  atom  recombination  coefficient  from  that  data  and  calculate  the  atom  hop  rate,  also. 
Defining  R  to  be  the  rate  of  conversion  from  ortho  to  para,  one  obtains: 


(16) 


R  = 


cFQ 
F  +  Q 


where  c  is  the  steady-state  atom  concentration,  F  is  the  atom  hop  rate  and  Q  is  the 
intrinsic  conversion  rate  for  a  fixed  atom  [13].  The  conversion  rate  has  a  sharp  maximum 
near  10  K,  3delding  a  peak  rate  Rp  =  8  x  10’'*  s*^  The  concentration,  c,  is  found  from  the 
equation 

(17)  c  =  N/Nt=  [K/2aN/]‘^' 
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By  dividing  R  in  Equation  16  by  the  peak  rate,  one  obtains  a  simple  function  of  the  ratio 
of  r/Q,  namely  R/Rp  =  2(r/Q)'^/[l  +  T/Q] .  Defining  =  T/Q,  the  value  of  x  can  be 
obtained  as  the  solution  of  a  quadratic, 

(18)  X  =  Rp/R  ±  V[(Rp/R)='  - 1] 

The  ratio  F/Q  can  be  obtained  at  each  temperature.  The  value  of  Q  is  obtained  from  the 
peak  rate  and  this  sets  the  scale  for  the  hop  rate  also.  The  values  of  F  obtained  by  this 
analysis  are  given  in  Table  2.  The  recombination  coefficient  is  then  calculated  from  F, 
using  Equation  10  in  the  form  2a  =  2voZ  F. 


Table  2.  The  Atom  Hop  Rate  (F),  Extracted  from  the  Ortho-para  Conversion  Rate  in  Tj, 
its  Reciprocal  (r),  and  the  Atom  Concentration  (c)  in  Parts  Per  Million  (ppm).  Tabulated 
at  Various  Temperatures 


T(K)  R(s-')  Rp/R  r/Q  r  (s'*)  t  (sec)  c,  (ppm) 


6 

l.llE-04 

7.2100E+00 

4.86E-03 

6.4 

1.23E-04 

6.5000E+00 

6.00E-03 

7.9 

1.72E-04 

4.6500E+00 

1.20E-02 

9.5 

2.87E-04 

2.7900E+00 

3.40E-02 

10.6 

5.42E-04 

1.4800E+00 

1.51E-01 

10.7 

5.62E-04 

1.4200E+00 

1.70E-01 

12.6 

3.32E-04 

2.4100E+00 

2.12E+01 

14 

1.12E-04 

7.1400E+00 

2.02E+02 

15 

8.32E-05 

9.6200E+00 

3.68E+02 

4.39E-02 

2.28E+01 

2.54E+03 

5.42E-02 

1.84E+01 

2.28E+03 

1.08E-01 

9.22E+00 

1.62E+03 

3.07E-01 

3.25E+00 

9.60E+02 

1.37E+00 

7.31E-01 

4.56E+02 

1.53E+00 

6.52E-01 

4.29E+02 

1.92E+02 

5.22E-03 

3.84E+01 

1.83E+03 

5.48E-04 

1.25E+01 

3.33E+03 

3.01E-04 

9.20E+00 

In  Figure  5,  the  recombination  coefficient  obtained  by  ESR  techmques  for  H  atoms 
in  crystalline  is  plotted  for  comparison  with  the  T  atom  recombination  coefficient 
obtained  by  modeling  the  NMR  data  on  solid  T2.  The  temperature  is  scaled  by  the 
quantum  parameter  for  H2,  namely  A  =  1.731  since  A  =  1  for  T2. 
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Figure  5 
a  Versus 

The  plot  symbols  symbols  are:  Tj  [Mn]  and  [v]. 


The  agreement  is  outstanding  even  though  the  recombination  coefficients  are  no  longer 
following  a  simple  thermally  activated  form. 

The  data  taken  in  Russia  on  atoms  in  and  Dj  [9,10]  represent  data  on 
amorphous,  not  crystalline  solids.  Up  to  now  there  has  been  an  attempt  to  understand 
the  two  sets  of  data  in  terms  of  the  extensive  database  that  exists  on  the  crystalline 
hydrogens.  A  different  approach  will  be  used  here.  Their  samples  were  prepared  by 
depositing  streams  of  or  Dj  gas,  containing  atoms  produced  by  microwave  discharge, 
onto  cold  substrates,  held  at  2  K  or  below.  The  recombination  coefficients  for  the 
amorphous  samples  are  compared  in  Figure  6  with  those  from  crystalline  samples.  Note 
that  there  is  not  good  agreement  even  between  the  amorphous  and  Dj  experiments, 
nor  is  there  agreement  between  the  amorphous  and  crystalline  samples.  There  remains 
the  possibility  that  the  samples  never  reached  equilibrium  and  the  vacancies  were  not 
thermal  vacancies. 
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Figure  6 
a  Versus 


2.2.  Suppression  of  Recombination  of  Atoms  in  3-D  Hydrogen  Hosts 

In  its  simplest  form,  the  diffusion  coefficient  D  can  be  written  as; 

(19)  D(T)  =  Dox, 

where  Dq,  at  low  temperatures,  would  be  the  diffusion  coefficient  appropriate  for 
activation  over  the  potential  barrier  between  molecules  and  is  the  concentration  of 
vacancies.  As  a  first  approximation,  it  is  found  that  x^  =  x^  -  (Vi)x3  where  x^  is  the  initial 
concentration  of  vacancies  at  a  given  temperature  and  x^  is  the  concentration  of  atoms. 
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Thus  D(T)  =  Do(x^  -  (V2)xJ  so  that  as  the  atom  concentration  increases,  the  diffusion 
coefficient  (and  hence  the  recombination  coefficient)  decreases,  making  larger  atom 
populations  possible.  Using  this  modified  value  for  the  diffusion  coefficient  leads  to  very 
large  differences  in  the  interpretation  of  experimental  results.  Specifically,  it  predicts  that 
(1)  the  atom  growth  at  a  given  temperature  can  have  the  shape  observed  in  electron 
paramagnetic  resonance  (EPR)  experiments  of  Sharnoff  and  Pound  [22]  and  later  by 
Collins  [23];  and  (2)  atom  populations  larger  that  the  above  calculated  steady-state  values 
can  be  obtained.  The  strange  atom  growth  feature  has  defied  explanation  for  30  years, 
until  now. 

The  theoretical  work  reported  here,  done  in  collaboration  with  Dr.  Peter  Fedders, 
Department  of  Physics,  Washington  University,  has  produced  a  breakthrough  in  the 
understanding  of  the  kinetics  of  the  atom  growth  in  solid  hydrogen  containing  tritium. 
The  tritium  undergoes  a  beta  decay  with  a  rate  that  is  almost  constant  over  the  life  of  an 
experiment  since  the  decay  half-life  is  12.3  years. 

The  electron  emerging  from  the  beta  decay  loses  its  energy  by  ionizing  and/or 
exciting  hydrogen  molecules  in  its  path,  producing  atoms  of  hydrogen.  If  the  same 
processes  occurred  in  the  solid  as  are  observed  in  the  gas,  each  beta  decay  would  produce 
about  1000  atoms.  These  atoms,  once  produced,  are  metastable  and  decay  (through 
recombination  into  molecules)  with  other  atoms.  Since  the  production  of  atoms  via  the 
tritium  beta  decay  is  constant  in  time,  recombination  limits  the  atom  population  and,  in 
the  simplest  of  cases,  a  steady-state  population  of  atoms  results.  This  simple  case  is  not 
always  observed  and,  recently,  insight  was  obtained  into  the  process  that  may  make  it 
possible  to  maintain  larger  populations  of  atoms  than  previously  thought. 

If  n  is  defined  to  be  the  atom  density  (number  of  atoms  per  cm®),  then  the  simplest 
rate  equation  describing  the  production  and  recombination  of  atoms  is: 

(20)  ^ 

at 

where  is  the  constant  rate  of  atom  production  in  atoms  per  second  per  cm®.  This  is  a 
known  quantity  and  can  be  calculated  from  the  beta-decay  rate:  =  1,14  x  10‘^  (xt) 
decays  per  cm®  per  second;  where  x-j-  is  the  tritium  concentration  in  the  sample;  and  the 
number  of  atoms  produced  per  beta  decay,  approximately  1000;  yielding  Kj  «  1.14  x  10^^ 
(Xy)  atoms/cm®/s.  The  recombination  coefficient,  written  as  2a,  is  another  matter.  It  has 
been  measured  directly  in  only  one  experiment  and  is  less  well  understood.  In  a  graph  of 
n  versus  time  (t),  starting  from  n  =  0  at  t  =  0,  the  atom  growth  curve,  the  above 
expression  predicts  that: 


(21) 


n  -  «otanh 


where  n^  is  the  steady-state  value  of  the  atom  density,  =  Ky2a  and  t  =  The 

predicted  growth  curve  is  smooth  beginning  at  zero  and  saturating  at  the  value  ng. 
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The  two  experiments  [22,23]  that  have  observed  growth  curves  that  are  very 
different  from  the  predicted  curve  both  found  a  distinct  plateau  in  the  curve  after  an 
initial  buildup  followed  by  a  second  region  of  atom  growth.  This  behavior  is  shown  in 
Figure  7. 


n(PPM) 
linewidth  (G) 


Time  (min.) 

Figure  7 

Atom  Density  as  a  Function  of  Time 


One  approximate  form  for  the  recombination  coefficient  is  2a  =  47rRJD(T),  where 
is  the  distance  (from  another  atom)  where  recombination  is  certain  (usually  taken  to  be 
the  lattice  spacing)  and  D(T)  is  the  diffusion  coefficient  for  an  atom.  Since  D(T)  is 
determined  by  the  vacancy  motion,  filling  of  the  vacancies  reduces  the  recombination 
coefficient  making  the  second  atom  buildup,  after  the  initial  plateau,  possible.  This 
model,  without  the  feature  of  vacancy  filling,  has  been  previously  used  to  describe  the 
crudest  features  of  atom  dynamics,  namely  constant  production  (KJ  and  temperature 
dependent  two-particle  recombination. 

This  simple  rate  expression  agrees  in  its  time  dependence,  when  =  0,  with  the 
only  direct  measurement  [9]  of  the  recombination  coefficient  (2a).  This  simple  expression 
predicts  that  the  atom  density  increases  as  the  hyperbolic  tangent  and  saturates.  It 
cannot  yield  a  time  dependence  such  as  that  seen  by  Sharnoff  and  Pound  [22]  and  also  by 
Collins  [23]  under  a  variety  of  conditions.  Such  a  time  dependence  is  characterized  by  an 
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inflection  point  in  the  curve  of  n(t)  versus  t.  Thus  the  second  derivative  of  n(t)  versus  t 
must  be  equal  to  zero  at  an  inflection  point. 

Making  the  simple  modification  of  equating  the  recombination  coefficient  (2a)  to 
(2a)oX^,  where  is  the  concentration  of  vacancies,  produces  a  rate  equation  that  has  an 
inflection  point  and  can,  at  least  crudely,  describe  the  observed  atom  time  dependence. 
The  concentration  of  vacancies  is  roughly  given  by:  =  x^-  y.J2,  where  x^  is  the  atom 
concentration  given  by  n/N  with  N  the  number  density  of  molecules  and  x^  the 
equilibrium  concentration  of  vacancies.  Introducing  this  into  the  initial  rate  equation 
gives: 


(22) 


dt 


=  K,  -  (2a), 


The  second  derivative  is  zero  when  x^  =  (4/3)x^.  The  steady-state  solution  of  this  cubic 
equation  is: 


(23) 

where  -  Ky(2a)oX^,  the  atom  density  that  would  be  obtained  in  the  steady-state  if 
atom  production  didn't  fill  the  vacancies  and  inhibit  diffusion  and  recombination. 

If  X  =  Xj/X^  and  Xq  =  ny^x^  are  defined,  then  the  above  condition  for  a  steady- 
state  solution  becomes: 


The  right  hand  side  is  always  positive  and  the  left  hand  side  can  be  positive  on  the 
interval  x  =  0  to  x  =  2,  the  inflection  point  being  at  x  =  4/3.  The  maximum  difference 
between  the  two  terms  on  the  left  hand  side  is  16/27,  also  occurring  at  x  =  4/3.  Thus 
depending  on  the  ratio  of  the  saturation  concentration  of  atoms  to  the  concentration  of 
vacancies,  there  can  be  a  steady-state  solution  for  x  <  4/3.  For  such  a  solution,  the  atom 
concentration  will  be  approximately  equal  to  (or  less  than)  the  vacancy  concentration. 

There  is  one  important  implication  of  this  simple  result:  the  larger  the  number  of 
vacancies,  the  higher  the  atom  concentration,  in  the  steady-state  or  at  the  beginning  of  a 
non-equilibrium  buildup  where  there  is  an  inflection  point.  This  point  is  new  and  has 
strong  implications  for  preparing  hydrogen  solid  to  hold  the  maximum  number  of  atoms. 

When  the  quantity  Xq  is  very  small,  a  steady-state  solution  is  possible  with  x  =  Xq. 
The  quantity  Xq  is  very  temperature  dependent,  as  seen  from  the  expression: 
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(25) 


thus  at  higher  temperatures,  a  steady-state  solution  could  exist.  At  lower  temperatures 
could  exceed  its  critical  value  of  16/27  and  no  steady-state  solution  would  exist,  only  an 
inflection  point  and  a  monotonically  increasing  atom  concentration. 

By  defining  a  dimensionless  time,  t^  =  [N(2a)oXvo^]t,  the  rate  equation  becomes: 
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which  is  integrable.  The  time  dependence  predicted  by  this  equation  is  shown  in  Figure  8. 
Note,  when  Xq  exceeds  the  critical  value,  the  time  dependence  predicted  is  that  seen  by 
Shamoff  and  Pound  [22].  In  this  figure,  the  value  of  Xq^  is  taken  to  be  (16.1/27),  just 
above  the  critical  value  of  (16/27).  When  the  term  that  takes  vacancy  filling  into  account 
is  omitted,  the  scaled  concentration  (x)  increases  uniformly  and  then  saturates  at  about  x 
=  0.77.  When  vacancy  filling  is  taken  into  account,  the  scaled  concentration  increases 
uniformly  to  a  much  larger  plateau  value  before  a  runaway  region  of  increase. 
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Figure  8 

The  Predicted  Concentration  versus  Dimensionless  Time 


The  right-hand  side  of  Equation  26  can  be  defined  to  be  B(x).  B(x)  has  a  minimum 
at  X  =  4/3,  independent  of  Xq,  as  can  be  seen  in  Figure  9.  When  =  (16/27),  the 
minimum  value  is  exactly  zero  so  this  value  is  the  critical  value.  The  points  where  B(x)  = 
0  are  called  fixed  points.  They  become  singularities  of  the  following  integral  and  are  the 
limiting  values  of  the  solutions  to  the  differential  equation  for  the  atom  concentration. 


(27) 


When  Xq^  <  (16/27),  the  solution,  starting  at  x  =  0,  proceeds  smoothly  to  the  first  turning 
point,  the  smallest  root,  x^,  on  the  positive  x  axis  of  the  cubic  equation  0  =  Xg^  -  x^  +  0.5 
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X®.  In  simple  terms,  the  atom  concentration  would  build  smoothly  and  then  saturate  at 
the  value  Xj.  It  is  interesting  to  note  that  the  value  of  Xj  depends  on  Xq.  When  Xq^  = 
(16/27),  the  value  of  x  saturates  at  x  =  4/3.  This  model  is  capable  of  explaining  the  time 
dependence  seen  by  Sharnoff  and  Pound  [22]  and  predicts  that  more  metastable  particles 
can  be  trapped  in  a  solid  by  filling  the  vacancies  and  reducing  the  diffusion. 
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Figure  9 
B(x)  Versus  x 


2.3  Production  of  Atoms  In  Situ 

Atoms,  produced  and  trapped  in  solid  hydrogen  hosts,  have  been  studied  by  a 
variety  of  experimental  techniques  but  detailed  knowledge  is  still  quite  fragmented. 
Admittedly,  the  experiments  have  all  been  complicated,  but  just  a  few  clear  results  have 
been  established  for  a  system  that  has  the  potential  for  storing  large  quantities  of  energy 
(in  the  form  of  atoms). 

2.3.1.  Motivation  for  a  New  Model  of  Atom  Production.  Leach  and  Fitzsimmons 
[9]  produced  atoms  of  hydrogen  in  solid  molecular  hydrogen  by  bombarding  the  solid 
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with  160  keV  pulsed  electron  beams.  The  electron  pulses  produced  atoms  in  the  solid  that 
were  then  studied  from  the  time  decay  of  their  distinctive  ESR  signal.  From  their 
estimates,  300  atoms  were  produced  per  160  keV  electron.  Since  the  signal  intensity  is 
proportional  to  the  number  of  atoms,  they  were  able  to  conclude  that  the  atom  density 
decay  followed: 


(28) 
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From  which  one  deduces  a  rate  equation  for  the  atom  density. 
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previously  used  and  known  as  the  bimolecular  recombination  equation,  the  one  expected 
from  a  process  involving  two  entities.  In  addition,  they  found  that  the  recombination 
coefficient,  2a,  was  strongly  temperature  dependent,  with  an  activation  energy  (195  K) 
that  matched  that  of  the  self-diffusion  coefficient  of  molecules  in  solid  hydrogen.  The 
lower  the  host  temperature,  the  smaller  the  diffusion  coefficient  and  hence  the  smaller  the 
atomic  recombination  coefficient. 

Prior  to  the  Leach  and  Fitzsimmons  work,  Sharnoff  &  Pound  [22]  had  made  a 
careful  study  of  solid  deuterium  that  contained  nominally  1%  of  tritium  (Tj),  again  using 
ESR  techniques  to  detect  the  D  atoms  produced  by  the  tritium  beta  decay.  Their  atom 
density  as  obtained  from  their  signal  intensity  grew,  then  saturated,  then  grew  again,  all 
at  constant  temperature.  One  would  expect  that  atoms,  continuously  produced  by  the 
beta  decays,  would  eventually  reach  a  maximum  population  where  additional  growth 
tendencies  would  be  balanced  by  recombination.  The  atom  growth  curve  observed  by 
Sharnoff  and  Pound  has  never  been  satisfactorily  explained  until  this  study.  From  Leach's 
work,  there  should  be  temperature  dependence  in  both  the  growth  curve  and  the  final 
atom  density  because  of  temperature  dependence  of  the  recombination  coefficient. 
Between  4  K  and  1  K,  Sharnoff  and  Pound  did  not  observe  any  temperature  dependence. 
Leach  estimated  that  23  atoms  were  produced  per  tritium  beta  decay  in  the  Sharnoff  and 
Pound  experiment. 

Webeler  [18]  studied  containing  small  amounts  of  Tj  but  he  didn't  use  ESR 
techniques.  He  studied  the  heat  released  when  atoms  recombined  in  the  host  solid  and 
attempted  to  quantify  storage  times  for  the  energy  resident  in  the  trapped  atoms,  again 
produced  by  the  tritium  beta  decay.  Rosen  [16]  and  Zeleznik  [17]  provided  a  more 
detailed  analysis  of  the  Webeler  experiments.  Rosen's  approach  is  examined  here. 

In  his  most  general  approach,  Rosen  treated  two  types  of  atoms,  mobile  atoms, 
that  were  really  ballistic  and  trapped  atoms.  Mobile  atoms  were  produced  by  the  beta 
decay  of  tritium  with  about  3  eV  of  kinetic  energy  and  percolated  through  the  lattice  like 
thermal  neutrons,  eventually  becoming  trapped  after  traveling  on  the  order  of  1  mm.  The 
trapping  of  the  mobile  atoms  in  the  solid  led  to  the  second  type  of  atom.  While  either 
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type  of  atom  could  contribute  to  the  stored  energy,  the  trapped  atoms  outnumbered  the 
mobile  atoms  by  a  factor  of  based  on  Rosen's  choice  of  parameters.  Diffusion  was 
considered  but  then  discarded  because  it  complicated  the  analysis  and  didn't  seem  to  be 
required  to  explain  the  experiments.  Estimates  were  made  of  all  the  pertinent  parameters 
and  recombination  was  treated  by  gas-phase  calculations  of  the  cross-section  leading  to  a 
temperature-independent  recombination  coefficient  (2a  =  2.6  x  cmVs)  applicable  to 
mobile  atom  —  mobile  atom  events  or  mobile  atom  —  trapped  atom  events.  If  an 
oversimplified  version  of  Rosen's  model  is  used,  reduced  to  one  type  of  atom,  say  the 
trapped  atoms,  the  atom  density  would  obey  a  rate  equation  given  by  Equation  20,  where 
K,  would  represent  the  constant  atom  production  from  tritium  beta  decays  and  would  be 
roughly  2a/10^  Rosen  assumes  that  the  atom  density  is  uniform  across  the  sample 
which,  in  turn,  predicts  a  steady-state  atom  population  given  by  rio^  =  KJa.  This  buildup, 
starting  from  n  =  0,  and  saturation  is  shown  in  Figure  8,  the  curve  without  vacancy 
filling. 

Although  Rosen's  model  is  valuable  in  understanding  the  heat  released  in  Webeler's 
experiments  and  points  out  the  critical  role  of  the  coupling  between  the  sample  and  the 
refrigerator,  there  are  several  difficulties  with  it: 

(1)  the  model  vastly  overestimates  the  volume  within  which  the  atoms  are 
produced,  because  an  H  atom  does  not  travel  through  solid  hydrogen  like  a  thermal 
neutron  but  thermalizes  in  a  much  shorter  distance; 

(2)  there  appears  to  be  no  mechanism  for  untrapping  an  atom  since  3  eV  would 

be  required  to  convert  a  trapped  atom  into  a  mobile  one;  and 

(3)  the  model  cannot  explain  the  observed  temperature  dependence  of  the  many 
experiments  on  atoms  in  solid  hydrogen,  other  than  the  Webeler  experiment. 

First,  the  recombination  coefficient  has  been  directly  measured  by  Leach  [9]  in  an 
experiment  where  the  H  atoms  were  produced  by  an  electron  gun  and  the  ESR 
measurements  were  made  with  the  beam  turned  off.  He  found  that  a  =  6.0  x  10  cm  /s 
at  6.75  K  and  rapidly  decreased  with  temperature.  Russian  data  [10,11],  where  the  H 
atoms  were  produced  by  trapping  atoms  formed  in  a  microwave  discharge,  interpreted  by 
Souers  [24],  would  lead  to  values  of  a  =  1.5  x  10'^^  cmVs  at  2  K  and  a  =  7.5  x  10'^  cmVs 
at  1.4  K.  It  can  be  argued  that  both  of  these  experiments  are  sensitive  to  the 
recombination  of  trapped  atoms,  not  considered  by  Rosen.  However,  in  the  temperature 
range  of  interest  in  the  Webeler  experiments,  the  recombination  used  by  Rosen  differs 
from  the  estimates  made  in  the  solid  state  by  at  least  5  orders  of  magnitude.  For  a 
recombination  coefficient  of  the  magnitude  used  by  Rosen,  the  calculated  atom  density 
from  the  heat  observed  in  a  triggered  energy  release  could  not  be  established  for  the  atom 
production  term  he  used. 

Second,  Rosen  estimated  the  volume  for  trapping  of  atoms  and  multiplied  that 
volume  by  the  beta-decay  rate  to  compare  with  the  actual  volume  of  the  sample.  Since 
his  calculated  volume  per  second  receiving  new  mobile  H  atoms  was  comparable  to  the 
sample  volume,  he  concluded  that  it  was  appropriate  to  consider  the  volumetric  rate  of 
production  of  mobile  H  to  be  approximately  constant  and  uniform  throughout  the  Hj 
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sample.  In  the  new  model,  where  the  parameters  are  obtained  from  other  experiments 
where  possible,  a  much  smaller  critical  volume  (not  a  trapping  volume )  is  obtained  that 
leads  to  a  small  fraction  of  the  sample  receiving  new  atoms  per  second  for  the  beta-decay 
rate  appropriate  to  pure  Tj  solid.  In  this  case,  it  is  not  appropriate  to  consider  the  atom 
production  to  be  uniform  and  atom  production  will  be  treated  differently. 

A  new  model  is  proposed  that  is  simpler  in  some  respects  but  more  complicated  in 
others,  than  Rosen's.  This  model  attempts  to  explain  the  ESR  experiments  (that  Rosen 
was  not  addressing)  and  draws  on  other  experiments  such  as  the  measurements  of  the 
range  of  beta  particles  in  solid  hydrogen  for  its  parameters. 

The  pertinent  features  of  this  new  model  are:  (1)  a  new  approach  to  atom 
production,  where  the  solid  is  treated  as  consisting  of  regions  with  different  atom 
densities,  not  a  homogeneous  medium,  leading  to  a  different  time  dependence  of  the  atom 
density  at  fixed  temperature;  (2)  a  new  interpretation  of  the  ESR  results  is  made  possible; 
(3)  the  model  has  some  predictive  ability  and  can  give  an  upper  bound  for  the  maximum 
atom  density  for  a  specified  production  process;  (4)  it  predicts  that  the  production  term, 
used  in  the  model,  can  depend  on  the  materials  properties  of  the  host  and  possibly  the 
temperature  and  not  just  the  fraction  of  tritium  in  the  sample;  and  (5)  it  can  guide 
experiments  attempting  to  make  large  atom  densities  by  electron  bombardment. 

In  the  analysis  of  most  of  the  experiments  to  date,  the  level  of  detail  of  the 
observations  has  not  warranted  using  a  theoretical  picture  with  the  (unphysical) 
complexity  of  Rosen's  (i.e.,  two  types  of  atoms)  so  previously  the  problem  was  treated  as 
though  only  one  type  of  atom  existed.  In  trying  to  reconcile  the  atom  densities  obtained 
from  a  model  analysis  of  ortho-para  conversion  times  in  Tj  (or  DT)  [25]  with  (1) 
measured  atom  densities  in  ESR  experiments  [26]  and  (2)  the  time  dependence  of  the 
density  after  a  temperature  change  (or  heat  spike)  [27],  it  is  apparent  that  the  single  t3^e 
of  atom  model,  with  homogeneous  production  of  atoms,  that  has  been  used  is  inadequate 
to  explain  the  data. 

The  beta  decay  of  a  single  tritium  atom  creates  atoms  in  a  restricted  volume  (a 
bubble)  of  the  crystal.  In  this  region,  there  can  be  a  dramatic  thermal  response  (local 
heating)  and  atom  recombination.  Atoms  created  in  this  bubble  can  recombine  there  or 
diffuse,  enlarging  the  bubble.  The  critical  initial  bubble  volume  is  estimated  from  the 
measured  range  of  electrons  in  solid  and  a  calculation  of  the  distance  traveled  by  a 
20  eV  secondary  electron,  making  the  bubble  shape  that  of  a  long,  thin  cylinder.  An 
effective  atom  production  rate  is  obtained  that  depends  on  the  temperature  and  the 
material.  Application  of  this  model  to  electron  paramagnetic  resonance  experiments 
indicates  that  many  atoms  can  go  undetected  in  ESR  measurements  when  the  atom 
density  is  large.  The  bubble  model  is  used  to  estimate  the  low  temperature  ortho  to  para 
conversion  rate  in  solid  tritium. 

2.3.2.  A  New  Model  of  Atom  Production  in  Solid  Hydrogen  Hosts.  Atoms  can  be 
produced  in  situ  by  electron  bombardment  of  Hj  or  by  radioactive  decay  (beta  decay)  of 
tritium  nuclei  incorporated  into  the  host  solid  (or  even  forming  the  host  solid).  In  either 
case,  the  role  of  an  electron  is  crucial.  The  motion  of  an  electron  injected  into  solid  Hj  is 
complicated  and  can  only  be  described  in  a  statistical  sense.  Historically,  this  problem  has 
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been  the  subject  of  interest  to  high  energy  physics  but  for  electron  energies  much  above 
those  of  interest  for  tritium  beta  decay  or  for  the  production  of  atoms  in  solid  Hj. 

A  collision  with  a  nucleus  would  lead  to  large  angular  deflection  while  a  collision 
with  another  electron  bound  to  a  molecule  could  result  in  almost  no  angular  deflection. 
The  former  situation  leads  to  multiple  Coulomb  scattering  and  for  electrons  with  energies 
of  400  keV,  incident  on  solid  argon,  10%  of  the  energy  is  deposited  behind  the  starting 
plane.  For  electrons  incident  on  solid  Hj,  by  the  time  the  energy  has  fallen  to  1.4  MeV, 
there  has  been  a  one  radian  deflection  but  this  effect  has  not  been  quantified  in  Hj.  The 
latter  situation  results  in  an  accumulation,  as  the  electron  loses  its  energy,  of  a  large 
number  of  small  angular  displacements.  These,  also,  have  not  been  quantified  in  solid  Hj. 

The  approach  followed  differed  from  that  of  Rosen  [16]  in  that  the  focus  was  on  a 
single  beta  decay  and  its  effect  on  the  surrounding  region  of  the  solid.  The  region  created 
by  the  beta  decay  was  labeled  a  bubble  and  the  sudden  increase  of  energy  available  in  this 
region  may  form  a  hot  spot.  A  more  careful  treatment  would  take  into  account  the  actual 
spectrum  of  tritium  beta-decay  energies  but  only  the  decays  characterized  by  the  mean 
energy  of  5.7  keV  were  considered. 

The  energy  loss  of  an  electron  in  matter  has  been  successfully  described  by  a 
continuous  slowing  down  approximation  where  it  is  assumed  that  the  electron  has  many 
collisions  (380  for  an  initial  energy  of  5.7  keV)  with  the  atomic  electrons  and  multiple 
Coulomb  scattering  is  ignored.  Although  the  electron  does  not  suffer  that  many  collisions 
and  it  does  not  transfer  its  energy  in  very  small  (but  nearly  identical,  15  eV)  bits,  the 
predictions  that  emerge  are  consistent  with  the  limited  experimental  observations.  The 
prediction  of  the  electron  practical  range  in  and  its  dependence  on  the  initial  electron 
kinetic  energy  is  a  good  example. 

For  this  initial  energy,  the  trajectory  of  the  electron  is  approximately  straight  but 
many  secondary  electrons  are  produced  as  the  initial  electron  slows  down.  Accordingly, 
the  region  influenced  by  this  electron  (the  bubble)is  presumed  to  be  a  cylinder.  One 
dimension  of  the  bubble  volume  was  chosen  to  be  the  practical  range  (L)  of  a  5.7  keV 
electron  in  solid  hydrogen  and  the  transverse  dimension  (R)  was  taken  to  be  the  range  of  a 
20  eV  secondary  electron  produced  in  large  numbers  from  the  energy  loss  of  the  primary 
beta  electron. 

Schou  and  Sorensen  [28]  measured  the  range  (L)  of  electrons  in  H2  and  D2  for  initial 
kinetic  energies  (E)  between  500  eV  and  3000  eV  by  using  thin  films  of  hydrogen.  They 
obtained  the  result  that: 

(30)  L  = 

where  K  is  a  constant  and  p  is  the  density  (of  U2,  D2,  or  T2).  Thus  an  electron  with  initial 
energy  of  5.7  keV  would  have  a  range  in  solid  H2  of  3.9  //m  or  in  solid  T2a  range  of  3.2  ^m 
before  losing  all  its  energy.  According  to  this  picture,  the  electron  would  lose  its  energy  in 
15  eV  increments  and  have  about  380  collisions  before  stopping.  In  fact,  the  number  of 
actual  collisions  could  be  much  smaller. 

The  transverse  dimension,  the  range  (RJ,  can  be  estimated  from  the  known 
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stopping  power  formula.  From  Souers  [24],  (dE/dx)  =  5.21  x  10  ‘°  J/m,  for  20  eV  initial 
energy  and  with  the  choice  of  5  eV  as  the  typical  (but  arbitrary)  energy  loss  per  collision 
in  the  low  energy  range.  Thus: 
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for  and  12.5  A  for  Tj  (because  of  the  larger  density).  Doubling  the  energy  loss  per 
collision  to  the  value  more  appropriate  for  dissociation  doubles  the  range  but  will  make  no 
appreciable  change  in  the  estimates  given  by  the  bubble  model,  since  the  effective  atom 
production  rate  depends  only  on  the  logarithm  of  the  transverse  dimension  so  that  a 
factor  of  3  error  in  would  produce  only  a  10%  change  in  the  rate. 

Thus  with  these  dimensions,  the  bubble  is  a  long,  thin  cylinder.  The  volume  of  a 
bubble  in  Tj  with  this  radius  (12.5  A )  and  length  3.2  //m  would  be  1.57  x  10  *^  cm^  while 
the  comparable  volume  in  would  be  2.9  x  10'‘^  cm®.  Given  the  beta  decay  rate  of  = 
1.15  X  10‘^  events/cm®/s,  in  a  pure  Tj  sample  of  volume  1  cm®,  1.8  x  10  ®  cm®  of  new 
bubbles  are  created  each  second.  For  the  2%  Tj  samples,  such  as  those  studied  here,  there 
are  3.6  x  10'®  cm®  of  new  bubbles  formed  each  second  or  in  a  sample  of  volume  1  cm®,  only 
36  parts  per  million  (ppm)  of  the  sample  would  receive  new  atoms  each  second,  justifying 
use  of  a  model  based  on  inhomogeneous  atom  production. 

Since  the  bubble  volume  is  so  small,  the  thermal  response  could  be  unusual.  The 
average  beta  particle  leaves  5.7  keV  in  the  bubble  region  while  producing  about  800 
atomis,  if  the  gas  phase  result  is  used  for  atom  production.  As  soon  as  the  atoms  have 
thermalized,  they  store  400  times  4.5  eV  or  1.8  keV,  leaving  3.9  keV  to  be  dissipated  in  the 
bubble  and  its  vicinity. 

The  enthalpy  of  Tj  has  been  estimated  from  measurements  and  comparison  with 
known  values  for  Hj  and  Dj.  At  the  critical  point,  the  enthalpy  is  1500  J/mol  (or  15.6 
meV  per  molecule).  If  the  above  value  is  used  for  the  cylinder  volume,  there  are  5  x  10® 
molecules  with  the  average  molecule  having  15.6  meV  enthalpy  at  the  triple  point  for  a 
combined  total  of  311  eV.  Thus  if  3900  eV  is  deposited  into  a  region  containing  5  x  10® 
molecules  (or  7.8  meV  per  molecule),  the  region  will  not  reach  the  critical  point. 

However,  smaller  estimates  of  the  bubble  volume  are  possible;  for  instance,  one  could 
assume  the  bubble  volume  is  400  times  the  volume  of  a  sphere  of  radius  12.5  A.  This 
reduces  the  volume  by  a  factor  of  5  so  that  only  10®  molecules  are  in  the  volume.  In  this 
case,  there  is  enough  deposited  energy  to  reach  (and  exceed)  the  critical  point.  However, 
based  on  a  classical  analysis,  the  thermal  relaxation  of  such  a  region  would  be 
unbelievably  fast,  on  a  time  scale  of  10'*^  seconds,  shorter  than  an  optical  phonon 
vibrational  time.  The  new  region  might  have  very  different  characteristics  than  the 
original  solid;  it  might  not  have  the  same  long  range  order. 

The  experiments  to  date  have  been  sensitive  only  to  time  scales  of  seconds  to 
minutes.  In  a  fraction  of  a  second,  all  the  atoms  that  are  going  to  be  produced  because  of 
a  single  beta  decay  have  been  produced  and  many  have  had  the  opportunity  to  recombine. 
Of  those  atoms  remaining  in  the  bubble,  some  diffuse  away  from  the  point  where  they 
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were  created,  enlarging  the  bubble.  Subsequently,  they  can  be  detected  in  an  electron  spin 
resonance  experiment.  Those  atoms  that  are  counted  in  the  detection  process  have  their 
ESR  frequency  address  located  within  a  magnetic  field  window  (aH)  centered  on  the 
particular  hyperfine-split  line  in  question.  This  window  depends  on  the  signal-to-noise 
ratio  and  the  sensitivity  of  the  spectrometer.  All  the  atoms  created  in  the  bubble  do  not 
instantly  appear  in  the  ESR  signal;  if  they  are  located  within  two  lattice  sites  of  another 
atom,  the  ESR  frequencies  of  both  atoms  will  lie  outside  the  detection  window  and  will 
not  be  counted.  Atoms  produced  in  the  high  density  regions  of  the  bubble  must  then 
diffuse  through  space  into  regions  of  lower  atom  density,  decreasing  their  local  magnetic 
field  and  ma^ng  them  detectable. 

To  quantify  these  effects,  a  microscopic  equation  for  the  recombination  time  for 
the  atoms,  including  the  dependence  of  its  variables  on  the  material  properties  of  the  host, 
was  developed.  For  two  atoms  to  recombine,  they  first  have  to  hop  within  a  nearest 
neighbor  distance  of  each  other  and  then  they  must  recombine.  This  last  step  may  be 
decidedly  nontrivial  since  the  repulsive  core  for  the  H  to  potential  is  wider  than  that 
for  the  Hj  to  potential.  Thus  there  could  be  a  substantial  repulsion  between  atomic  H 
nearest  neighbors  in  the  host.  In  any  case,  the  average  recombination  time  T  for  atoms 
is: 


(32) 


T  = 


+  t. 


Zc 


where  is  a  numerical  constant  of  order  unity,  Z  is  the  number  of  nearest  neighbors,  and 
l/x,  is  the  recombination  rate  for  neighboring  H  atoms.  The  quantitylAh  is  the 
probability  per  unit  time  for  an  H  atom  to  hop  a  specific  neighboring  site  and  c  is  the 
concentration  of  atoms. 

In  the  bubble  of  H  atoms  being  considered  there  are  N  atoms  in  a  volume  V.  Thus 
c,  the  concentration  of  atoms,  is  given  by  the  equation 


(33) 


c  = 


V 


where  Vg  is  the  volume  of  a  cell  containing  one  Hj  molecule.  Thus  V/Vg  is  the  number  of 
unit  cells  or  molecules  in  the  bubble.  The  differential  equation  governing  the  decrease 
in  the  number  of  atoms  in  the  bubble  is 


(34) 


dN  _  N 

dt  T 


where  the  factor  2  obtains  because  each  recombination  destroys  two  atoms.  Further,  V , 
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which  appears  in  this  equation  though  the  factor  c  in  the  recombination  time  T,  depends 
on  time  because  of  diffusion.  This  time  dependence  can  be  well  approximated  by  the 
equation 


where  is  the  length  of  a  cylinder  and  R^,  is  the  cylinder  radius.  The  physics  of  this 
equation  is  just  that  for  a  diffusion  process  <x^>  =  2Dt. 

Combining  these  equations,  the  form  to  be  solved  can  be  written 


(36) 


J_  .^  =  _  2  ^0 

dt  ■  t,F(0 


where  r  is  given  by 


(37) 


r  = 


Vq  is  the  volume  per  particle,  and  V(t)  is  given  by  Equation  35. 

Note  that  the  recombination  coefficient  is  given  in  terms  of  the  parameters 
previously  introduced  as 


(38) 


rv, 


a  = 


(^jT,  +  Zt^)  T, 


Now  Equation  36  is  solved  with  the  boundary  condition  that  N  =  Nq  at  time  t  =  0, 
yielding  a  solution  that  can  be  written  as 


(39) 


JL  =  1 

’  [1  ♦  PWl 


where  h(y)  is  given  by  the  equation; 


(40) 


h(y)  = 


ln(l  -  e)  +  e 


The  parameter  y  =  sx  where  x  is  the  dimensionless  time; 
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and  s  =  2a/Ro  depends  on  the  shape.  The  parameter  e  in  Equation  40  is  defined  as  e  =  1  - 
VLo-  In  Equation  40,  h(y)  is  zero  for  y  =  0  and  unity  for  y  infinite. 

The  factor  fe  in  Equation  39  is  given  by: 


(42) 


^  [  4/WoVoV  ln(l  -  e)  f  € 


The  first  factor  in  the  expression  for  C»  depends  on  the  initial  concentration  of  atoms  in 
the  bubble  but  is  shape  independent,  while  the  second  does  depend  on  the  shape.  For  a 
long,  thin  cylinder,  e  «  1  and  the  expression  for  C>  becomes: 


(43) 


Besides  the  bubbles  that  are  nucleated  by  beta  decays,  one  will  have  a  rather 
uniform  background  density  made  up  of  bubbles  that  have  grown  very  large  from  the 
diffusion  (and  lost  atoms  due  to  recombination).  A  sensible  but  arbitrary  way  to  divide 
things  is  to  consider  the  atoms  in  a  bubble  to  belong  to  the  uniform  background  density 
when  the  atom  density  of  the  bubble  is  comparable  to  the  background  density.  In  this 
case,  the  evolution  of  atoms  in  the  background  density  is  described  by  Equation  20  but 
with  an  effective  production  constant. 

From  the  above  equations,  one  can  immediately  deduce  the  effective  number  of 
atoms  that  each  beta  decay  produces.  At  times  long  compared  to  the  time  needed  to 
diffuse  across  the  bubble  radius,  =  Ro^  th/4a^  the  number  density  of  atoms  in  the 
bubble  N(t)/No,  is  comparable  to  the  uniform  background  density.  Then  one  can  consider 
the  bubble  should  merge  into  the  uniform  background  and  its  remaining  atoms  will  be 
added  to  the  uniform  background.  Note  that  almost  all  of  the  recombination  happens 
within  several  diffusion  times,  Tj.  The  number  of  atoms  left  is  given  by  Equation  39  with 
h(y)  equal  to  unity.  Thus  one  obtains  the  equation 


(44) 


_1 _ 

1  +  P 


This  analysis  also  yields  an  effective  production  constant  for  Equation  10.  The  relevant 
equation  is 
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(45) 


K 


K 

1  -H  p 


This  is  the  quantity  that  belongs  in  Equation  22  since  only  a  fraction  of  the  atoms  survive 
long  enough  to  be  incorporated  into  the  uniform  background.  The  Hj  sample  can  be 
pictured  as  containing  a  rather  uniform  density  of  atoms  (the  background)  along  with 
some  high  density  bubbles.  As  the  number  of  atoms  in  a  bubble  decreases  and  (especially) 
as  the  volume  of  a  bubble  increases,  the  atom  density  of  the  bubble  becomes  comparable 
to  the  density  of  the  uniform  background.  Thus  the  bubbles  supply  the  uniform 
background  with  a  rate  constant  that  is  adjusted  because  many  of  the  atoms  don't  survive 
long  enough  to  be  incorporated  into  the  uniform  background.  In  this  case,  the  appropriate 
production  term  is  K^jf. 

The  dimensionless  quantity  is  the  key  parameter  in  the  above  expression.  Note 
that  it  is  proportional  to  the  number  of  initial  atoms,  as  one  might  have  expected,  and  it 
also  depends  on  the  bubble  shape.  Thus  for  example,  one  could  imagine  changing  Ng  or 
V(0)  by  slightly  doping  the  Hj  and/or  injecting  electrons  with  a  variety  of  incident 
energies. 

The  effective  production  constant  is  dependent  on  the  material  properties  of  the 
host  (i.e.  H,  D,  or  T)  and  not  just  on  the  fraction  of  radioactive  T  atoms.  Further,  this 
constant  can  be  temperature  dependent  if  r  is  less  than  unity. 

Two  limits  of  the  effective  production  term,  the  limit  obtained  when  Th  >  >  ^  t, 
and  the  opposite  limit  obtained  when  z  t,  >  >  x,,  are  considered.  (1)  When  the  hopping 
time  is  much  longer  than  the  recombination  time,  is  evaluated  by  taking  take  N„  = 

800,  the  gas  phase  value,  a  =  3.75  A  for  Hj,  =  4  mm  for  Hj,  and  Rg  =  12.5  A.  Thus 
from  Equation  43, 


(46a) 
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0.119  r  =  17.2 


making  Kjff «  K/20. 


(2)  When,  however,  the  recombination  time  dominates  the  problem. 


(46b) 


llrKa] 

f4l 

0 

In 

-  1 

=  0.119  r  =  1.43 

h 

^  ^TCZig  ^ 

Thus  R  could  be  much  less  than  unity  and  the  effective  production  term  would  then  be 
the  actual  atom  production  term.  This  possibility  introduces  an  important  temperature 
dependence  and  isotope  dependence  into  the  problem.  It  is  estimated  x^/x,  «V2  at  1.4  K, 
making  R  =  0.7  and  »  K/2. 

One  can  estimate  the  hopping  times  from  the  diffusion  data  and  the  recombination 
times  from  a  calculation  of  the  time  needed  to  flip  the  spin  of  an  atom  adjacent  to  another 
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atom,  a  Tj  calculation.  Two  atoms  with  parallel  spins  are  prohibited  from  recombining. 
While  there  may  be  difficulties  in  their  becoming  near  neighbors  in  the  lattice  due  to  an 
effective  repulsion,  once  they  are  near  neighbors,  they  could  recombine  almost 
instantaneously  on  the  timescale  of  the  ESR  experiments  if  their  spins  were  antiparallel. 
Thus  an  estimate  of  the  time  required  for  a  spin-flip  of  one  of  the  atoms  gives  insight  into 
the  time  constant  t,,  the  recombination  time. 

Abragam  [32]  has  estimated  the  transition  probability  per  unit  time,  P,  for  such  a 
spin-flip  using  a  Raman  process.  At  temperatures  T  <  <  0,  where  0  is  the  Debye 
temperature,  P  is  given  by: 


(47) 


P  «  732.5 


{  81 

f  rj^\ 

i  2  J 

[mv^l 
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where  Fj  is  proportional  to  the  second  derivative  of  the  dipole-dipole  energy,  m  is  the 
atom  mass,  Qq  is  the  Debye  frequency,  and  v  is  the  sound  velocity.  Rough  estimates  of 
this  rate  give  for  H  atoms,  P  =  (1.35  x  10’'*  s'*)T^  so  that  >  t,to  below  2  K. 

From  estimates  of  the  model  parameters,  it  is  found  that  the  initial  concentration 
of  atoms  in  a  bubble  is  about  5%.  Such  a  density  of  free  electronic  spins,  at  a  temperature 
of  about  4  K,  randomly  distributed,  would  represent  an  extremely  interesting  magnetic 
system.  There  are  no  existing  studies  of  the  magnetic  behavior  of  such  a  system  so  the 
arguments  below  will  be  based  on  extrapolations  of  results  from  other  experiments. 

What  are  the  implications  for  obtaining  even  greater  atom  densities^  In  one 
second,  1.15  x  10‘'‘new  bubbles  per  cm^  are  created,  with  each  bubble  initially  having 
volume  6  x  10‘‘’  cm^  Thus  a  small  fraction,  7  x  10  ®,  of  the  sample,  each  second,  receives 
new  atoms.  A  bubble  is  incorporated  into  the  uniform  background  in  a  time  of  order  the 
diffusion  time.  For  a  diffusion  time  of  150  s  (the  estimate  at  4  K  for  Tj),  only  1%  of  the 
sample  is  in  this  high  density  phase.  This  suggests  that  it  may  be  possible  to  achieve 
much  higher  atom  densities  in  the  solid  host,  but  not  higher,  for  the  calculation  based  on 
beta  decay  in  solid  Tj,  than  the  initial  atom  density  of  5  %. 

In  the  above  analysis,  it  was  assumed  that  a  single  electron  of  average  energy  of  5.7 
keV  can  create  1000  atoms  in  a  small  volume  (a  bubble)  estimated  to  be  of  order  6  x  10'‘® 
cm®.  The  bubble  expands  due  to  diffusion  and  atoms  in  the  bubble  recombine.  When  the 
atom  density  decreases  by  a  factor  of  10  from  its  initial  value,  only  half  of  the  atoms  have 
recombined,  the  major  effect  being  a  volume  increase  by  a  factor  of  5.  Thus  a  maximum 
atom  density  of  about  2%  does  not  seem  impossible. 

To  enhance  the  chances  of  obtaining  such  high  densities,  it  seems  that  limiting  the 
effects  of  recombination  by  increasing  the  applied  magnetic  field  and  possibly  increasing 
the  pressure  are  most  likely  to  succeed.  In  the  experiments  of  Smith  and  Squire  [33]  on 
solid  molecular  Hj,  modest  increases  in  the  pressure  (p  =  1  atm  to  p  =  230  atm)  made 
significant  increases  in  the  activation  energy  (E^  =  175  ±25  K  to  E,  =  280  ±25  K)  for 
diffusion  and,  by  implication,  the  diffusion  time.  Without  additional  data,  one  cannot 
speculate  on  the  modifications  required  for  this  bubble  model  due  to  changes  in  pressure 
and  magnetic  field. 
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2.4  Hj  Mcaiolayets  Absorbed  on  Activated  Carbon  Fibers 


When  the  spins  are  confined  to  lattice  positions,  the  NMR  lineshape  is  essentially 
Gaussian  in  nature  (Gaussian  FID)  and  its  width  (Aw)  is  well  characterized  by  the  second 
moment,  Mj,  with  (Aw)  »  / (Mj).  Time  dependence  (motion  or  spin  flips)  of  the  spins 
results  in  a  narrower  line  and  a  change  in  lineshape  from  Gaussian  to  Lorentzian 
(exponential  FID),  effects  predicted  by  Anderson  and  Weiss  [31].  If  the  motional 
frequency  is  F,  the  motionally  narrowed  linewidth  is  given  by:  (Aw)  «  a  result 

valid  so  long  as  F  >  V (Mj).  The  motion  responsible  for  this  narrowing  can  come  from 
actual  translational  hopping,  in  which  case  the  data  can  be  used  to  infer  diffusion 
coefficients,  or  it  can  be  caused  by  intramolecular  interactions  between  neighboring  oHj 
molecules  that  produce  spin-flips,  without  any  translational  motion.  These  latter 
interactions  are  absent  in  pure  HD  so  that  a  narrowed  HD  resonance  line  is  indicative  of 
translational  hopping. 

In  experiments  with  pure  HD  (about  1.9%  Hj  and  0.3%  Dj),  a  Lorentzian  lineshape 
was  observed,  indicative  of  motional  narrowing,  and  the  weakly  temperature-dependent 
observed  widths  (2.2  x  10^  at  13.6  K  to  4.4  x  10^  rad/s  at  5  K)  are  smaller  than  the  width 
predicted  from  the  second  moment  (5  x  10^  rad/s).  This  is  shown  in  Figure  10.  At  7.6  K, 
the  hopping  frequency  calculated  from  the  data  is  F  =  8  x  lO"  rad/s  and  the  2-D  diffusion 
coefficient  is  D  =  <a^>F/4  =  2.3  x  10‘”  cmVs,  comparable  to  the  3-D  value  for  Hj  at  11  K. 


A  Aw-Lorentz 

T  a-gaussian 
•  Aw  (from  M2) 


Figure  10 

linewidth  versus  T 
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Au)  (rad/s) 


The  NMR  lineshapes  of  Hj  molecules  absorbed  on  ACF  are  also  Lorentzian  and,  at 
the  starting  oHj  concentration  of  Xj  =  0.75,  all  widths  are  narrower  than  predicted  from 
the  calculated  value  of  Mj.  As  the  samples  age,  at  fixed  temperature,  their  linewidths 
increase,  exactly  the  opposite  behavior  expected  from  spins  on  a  rigid  lattice.  This  is 
shown  in  Figure  11. 


linewidth  Versus  Ortho  Concentration 
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This  would  indicate  that  the  hopping  frequency  decreases  with  decreasing  oHj 
concentration.  Nevertheless,  the  hopping  frequencies  are  comparable  to  those  found  in 
HD  and  the  diffusion  coefficients  are  of  order  of  magnitude  10'*^  cmVs.  The  weak 
temperature  dependence  of  the  hopping  frequencies  for  Xj  =  0.75  is  characteristic  of  an 
activation  energy  of  5  K. 

Since  the  NMR  signal  height  for  Hj  is  directly  proportional  to  the  number  of  oHj 
molecules,  the  decay  of  the  signal  at  fixed  temperature  can  be  used  to  study  the  0H2-PH2 
conversion  process.  Two  types  of  conversion  have  been  studied  previously:  intrinsic 
conversion  where  the  intermolecular  interaction  between  molecules  provides  the  magnetic 
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field  gradient  needed  for  conversion  and  catalyzed  conversion  where  a  paramagnetic 
impurity  such  as  an  H  (or  O)  atom  or  Oj  molecule  provides  the  needed  magnetic  field 
gradient. 

Using  a  strong  collision  approach  to  calculate  the  local  transition  probabilities,  one 
obtains  the  following  differential  equations  for  intrinsic  conversion: 

=  -  K  if  >  1 

T 

=  -  (N,)\Zv,W,X  if  <  1 


(48) 


dN, 

(a)  - i- 

^  ^  dt 

dN, 

(b)  - 1 

dt 


where  t(t  =  is  the  time  to  hop  to  any  of  the  Z  nearest  neighbor  sites,  Vq  is  the  volume 
per  particle  such  that  VqN-j-  =  1,  and  f^  is  a  correlation  function  that  gives  the  fraction  of 
near  neighbor  sites  at  the  new  site  that  were  not  near  neighbor  sites  at  the  old  site.  For 
the  hep  lattice,  Vq  =  aVV*2.  To  modify  this  for  2-D,  replace  Vq  by  A^,  the  area  per  particle, 
where  Aq  =  (/3/2)a^.  For  two  dimensions  f^  is  approximately  0.50. 

For  intrinsic  conversion,  (Nj)'^  is  linear  in  time  and  the  rate  has  been  measured  [32] 
for  Hj  and  found  to  be  1.9%  per  hour  or  5.3  x  lO'**  s'*.  A  careful  calculation  of  the  basic 
rate  in  has  been  done  [33]  finding  a  rate  of  1.67%  per  hour. 

Using  the  same  procedure  as  above,  the  differential  equations  for  catalyzed 
conversion  (with  c^  the  concentration  of  unpaired  electron  spins)  were  obtained: 


(«)^  =  -  (N,)  if  wfa)  T  >  1 

at  X 

(49) 

dN, 

(P)-^  =  -  (N,)  if  Wfa)  X  <  1 


In  both  cases,  there  is  an  exponential  decay  of  the  number  of  ortho  molecules,  Ni(t)  = 
Ni(0)  exp(-kt),  where  k  depends  on  which  of  the  above  limits  applies. 

By  using  the  factors  in  the  calculation  (but  using  the  experimental  value  of  the 
rate),  it  is  possible  to  predict:  the  intrinsic  conversion  rate  for  an  H2  monolayer  (2-D 
intrinsic  rate  Wjq  «  1.7  x  10'®  s'*)  and  the  transition  probability  per  unit  time  due  to  an 
unpaired  electron  spin  at  the  nearest  neighbor  spacing  (a).  For  3-D  Hj,  this  is  Wfa)  = 

5.46  X  10'^  s'*  (in  agreement  with  the  estimate  of  Shevtsov,  et  al.  [34]).  This  latter  rate  can 
be  adjusted  to  the  monolayer  spacing  giving  for  2-D,  ACF,  the  value  Wc(a)^^  »  0.177  s'*. 
The  measured  0H2-PH2  conversion  rates  are  about  a  factor  of  10  larger  than  the  predicted 
intrinsic  conversion  rates  and  the  signal  decay  is  approximately  exponential,  indicative  of 
catalyzed  conversion.  This  is  not  surprising  since  O2  contamination  of  the  surface  is 
assured  given  the  transport  and  handling  conditions  of  the  ACF.  Thus  intrinsic 


39 


conversion  was  ignored  in  the  analysis  here. 

There  is  an  additional  feature  to  the  time  decay  of  the  NMR  signal  at  fixed 
temperature,  namely,  the  conversion  rate  appears  to  increase  rather  abruptly  with  time,  a 
phenomenon  observed  previously  for  3-D  bulk  Hj  by  Schmidt  [32]  at  1.57  K  and  by  Rail, 
etal.  [35]  for  Hj  absorbed  on  zeolite.  This  is  shown  in  Figure  12. 
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Figure  12 

Ortho  Concentration  versus  Time 


In  the  fast  diffusion  limit,  the  rate  ki  is  given  in  Equation  49b:  kj  =  For 

Z  =  6,  fj  =  0.5,  and  W£(a)  =0.177  s'^  it  is  obtained  from  <ki>  that  c^  =  4  x  lO"*  (or  400 
ppm).  This  limit  implies  that  F  is  larger  than  0.177  rad/s.  Motionally  narrowed  NMR 
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lineshapes  require  Fs  larger  than  10'*  rad/s. 

In  the  slow  diffusion  regime,  then  kj  =  F  and  the  mean  rate  implies  that  the 
product  Ca  r  =  7x  10'® s'*,  is  not  known  independently  but  choosing  F  «10^  rad/s, 
consistent  with  the  NMR  line  narrowing,  would  imply  that  c^  «10'®,  an  unlikely  value 
(much  too  small  for  samples  exposed  to  air  in  shipping).  Thus  it  is  concluded  that  the 
data  is  in  the  fast  diffusion  limit  and  the  hopping  frequencies  can  be  determined  from  the 
NMR  line  narrowing  but  not  the  ortho-para  conversion  rates.  The  Hj  data  is  consistent 
with  the  HD  data  and  yields  diffusion  coefficients  on  the  order  of  D  «10'**  cmVs. 

There  is  no  explanation  for  the  increase  in  the  conversion  rate  observed  as  the 
samples  age  at  fixed  temperature.  Other  workers  [35,38]  attributed  the  increase  they 
observed  to  clustering  of  the  ortho  molecules  driven  by  the  electric  quadrupole-quadrupole 
interaction.  While  this  might  be  true  for  their  experiments,  it  is  not  a  likely  source  of 
increased  conversion  when  the  conversion  is  catal3/2ed  by  unpaired  electron  spins.  From 
the  NMR  lineshape  narrowing,  there  is  no  indication  of  increased  diffusion  as  the  samples 
age;  rather  the  opposite  is  true.  The  frequency  F  appears  to  decrease  (and  hence  D)  with 
decreasing  ortho  concentration  (Xj)  or  increasing  age. 
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3.0  SUGGESTIONS  FOR  FUTURE  WORK 


The  suggestions  for  future  work  are  divided  into  three  categories; 

(1)  conventional  diffusion  measurements  but  on  vapor  deposited  soHds,  (2)  non- 
conventional  diffusion  measurements  that  can  pinpoint  the  actual  molecule  hopping  time 
unambiguously  at  low  temperatures,  and  (3)  basic  measurements  of  the  spin-lattice 
relaxation  time  on  HD  and  molecules  in  monolayers. 

3.1  Conventional  Diffusion  Measurements 

Additional  measurements  of  the  self  diffusion  coefficient,  D,  in  vapor  deposited 
hydrogen  solids  (a-Hj),  from  near  the  triple  point  temperature  to  the  lowest  practical 
temperatures  using  conventional  NMR  techniques,  should  be  made.  There  is  now 
comprehensive  literature  on  these  same  measurements  in  solid  prepared  by 
crystallization  of  the  liquid  (c-Hj).  Thus  a  comparison  of  the  two  sets  of  data  will  make  it 
possible  to  quantify  the  differences  in  vacancy  concentrations  in  the  two  types  of 
samples. 

3.2  Diffuaon  Measurements  in  the  Rotating  Frame 

The  direct  coimection  between  molecular  motions  and  the  linewidth  and 
relaxation  of  magnetic  S5^tems  provides  one  of  the  most  powerful  techmques  known  for 
studying  molecular  motion.  Molecular  motion,  of  the  diffusive  type,  is  characterized  by  a 
mean  time  a  molecule  resides  at  a  given  site  between  jumps,  t  =  1/T.  For  the 
conventional  motional  narrowing  approach,  one  can  measure  t  as  long  as  t  is  less  than 
(T2)rl-  For  t5q)ical  hydrogen  solids,  (Tj)^.  is  determined  by  intermolecular  dipolar 
interactions  and  is  on  the  order  of  magnitude  of  30  to  400 

The  alternate  technique  first  discovered  in  1964  by  Slichter  and  co-workers  [36, 37] 
extends  the  range  of  measurement  of  motional  correlation  times  by  many  orders  of 
magnitude,  making  molecular  motion  quantifiable  when  x  <  Tj,  the  spin-lattice 
relaxation  time.  By  this  extension,  it  is  possible  to  observe  hopping  times  as  slow  as  one 
lattice  site  hop  per  hour.  This  extends  the  range  of  t  measurements  from  10"*  s  to  10^  or 
10®  s  (or  more),  a  factor  of  more  than  10*^.  This  is  just  the  extension  needed  to  study  ultra- 
slow  molecular  hops  in  solid  in  the  temperature  range  of  interest  to  the  HEDM 
program. 

Slichter,  using  spin-locking,  performed  an  adiabatic  demagnetization  in  the 
rotating  frame  of  reference  where  typical  NMR  experiments  are  treated.  From  his  analysis 
and  experiments,  he  concluded  that  very  slow  motions  could  be  detected  when  the 
motional  correlation  time  x„  was  smaller  than  Tj. 

This  technique  was  applied  to  NMR  experiments  in  by  the  Gaines  group  [38] 
and  in  Dj  by  Meyer  and  co-workers  [6] .  In  an  elegant  set  of  experiments,  they  were  able 
to  measure,  by  conventional  means,  the  motional  narrowing  approach,  x's  that  ranged 
from  10’**  s  to  2  X 10"*  s.  However  by  measuring  the  relaxation  time  in  the  rotating  frame 
of  reference,  they  extended  the  measurements  of  x  from  1  ms  to  10  s  at  temperatures  as 
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low  as  8.9  K  in  Dj.  By  a  scaling  argument,  this  would  correspond  to  a  temperature  of  5.8 
K  in  Hj.  Since  this  new  technique  permits  measurement  of  t  if  x  <  Tj,  it  is  important  to 
work  in  solid  hydrogen  where  the  relaxation  time  Tj  is  long. 

It  is  possible  however  to  obtain,  without  great  difficulty,  long  relaxation  times  in 
solid  HD  since  the  HD  molecules  have  no  intrinsic  relaxation  mechanism  of  their  own  — 
they  can  only  change  their  spin  orientation  if  they  have  oHj  molecules  in  their 
neighborhood,  making  cross-relaxation  possible.  There  it  is  necessary  to  either  purify  the 
HD  to  remove  oHj  impurities  (HD  samples  have  been  made  here  by  double  distillation  to 
remove  Hj  molecules)  or  to  hold  the  HD  sample  at  low  temperatures  for  a  few  weeks  and 
allow  the  oHj  impurities  to  convert  to  pHj,  which  pose  no  problem.  In  either  case,  it  is 
possible  to  produce  HD  samples  with  TiS  of  several  hundred  seconds. 

Because  HD  molecules  have  a  different  quantum  parameter,  the  temperature  of  5  K 
in  solid  Hj  corresponds  to  a  temperature  of  (5  K)  (250/197)  =  5.3  K  in  HD.  Thus 
experiments  on  HD  at  6.3  K  will  provide  the  needed  information  about  the  hopping  times 
in  Hj  down  to  5  K. 

3B  Measurements  of  the  NMR  Relaxation  Times  in  Hj  Monolayers 

Both  Hj  and  HD  molecules  confined  to  a  two-dimensional  array,  such  as  that 
found  for  monolayers  on  ACF,  are  in  mixed  rotational  states.  The  HD  molecules  are  no 
longer  in  a  pure  J  =  0  ground  state  but  have  a  significant  J  =  1  component;  even  the  p-Hj 
molecules  are  in  a  state  that  mixes  the  rotational  states  J  =  0  and  J  =  2.  This  miying  of 
the  states  nullifies  one  of  the  major  simplifications  found  in  3-D  solid  H2,  namely  that 
only  the  oHj  molecules  (J  =  1)  possess  a  non-zero  average  value  for  the  electric 
quadrupole-quadrupole  moment. 

Because  of  the  mixed  states  in  2-D  Hj,  there  must  be  measurements  of  the  spin- 
lattice  relaxation  times  in  both  HD  and  Hj  to  understand  the  basic  interactions  between 
molecules  on  these  surfaces.  It  is  entirely  possible  that  the  motional  narrowing 
interpretation  used  to  deduce  the  diffusion  coefficients  of  HD  and  Hj  molecules  in 
monolayers  on  ACF  is  incorrect  and  that  the  observed  effect  is  due  to  the  electric 
quadrupole-quadrupole  interactions  between  the  molecules  and  not  to  molecular  hopping. 
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4.0  CONCLUSIONS 


New  measurements  of  the  molecular  diffusion  coefficient  in  HD,  D-T,  and  Tj  solids 
were  provided  and  the  new  data  was  interpreted,  along  with  previous  NMR  and  EPR  data 
to  extract  both  molecule  and  atom  hop  frequencies  in  the  hydrogen  isotopes.  The  vacancy 
exchange  mechanism  is  responsible  for  both  self  diffusion  of  molecules  and  that  of  atoms. 
Based  on  this  data  and  analysis,  the  following  specific  conclusions  can  be  made. 

Diffusion  in  these  crystalline  solids  (c-Hj,  etc.)  is  controlled  by  the  number  of 
vacancies  in  the  lattice.  This  is  in  agreement  with  the  theory  of  diffusion  in  hydrogen  by 
Ebner  and  Sung  [8]  .  The  non-radioactive  isotopes,  Hj,  HD,  and  Dj,  have  values  of  Dq  that 
are  very  regular  and  predictable  from  the  relationship  Dq  =  (0.0241)A'^^  The  radioactive 
isotopes  have  much  larger  values  of  Dq  that  compensate,  in  part,  for  their  much  larger 
activation  energies.  The  apparent  activation  energies  may  arise  from  lattice  damage 
caused  by  the  radioactivity  that  acts  to  suppress  hopping  and  diffusion. 

Values  of  the  vacancy  formation  energy,  E^,  the  barrier  height  energy,  Et,  and  the 
energy  of  the  first  tunneling  level  in  the  hydrogen  potential,  £„  defined  by  Ebner  and 
Sung,  have  been  estimated  for  all  the  isotopes.  The  activation  energy  according  to  the 
Ebner  and  Sung  theory  is  given  by:  E3  =  E^  +  E,,  -  Ej.  Model  dependent  values  for  these 
quantities  have  been  given. 

The  values  for  the  microscopic  hopping  rates  of  the  molecules  were  obtained  from 
the  data.  The  vacancy  hopping  rate,  at  the  triple  point,  for  all  the  isotopes  was  calculated 
and  found  to  be  essentially  the  same  for  all  isotopes.  The  microscopic  atom  hopping  rates 
for  crystalline  samples  are  slightly  larger,  but  nearly  identical,  to  the  corresponding 
hopping  rates  for  molecules.  Thus  by  studying  the  molecule  hopping  rates  (say  via  NMR 
measurements)  one  can  infer  the  atom  hopping  rates. 

The  recombination  coefficients  (and  hence  hopping  rates)  for  crystalline  solids 
differ  from  those  of  amorphous  solids.  Extrapolation  of  the  data  on  crystalline  solids  to 
higher  temperatures  shows  that  atom  diffusion  is  also  dominated  by  vacancies  and  that  at 
the  triple  point  the  atom  hop  frequency  is  somewhat  larger  than  the  molecule  hop 
frequency. 

One  method  of  recombination  suppression  explains  the  bizarre  features  of  the 
buildup  curve,  possibly  pointing  the  way  to  storing  much  larger  densities  of  atoms  in  solid 
hydrogen.  It  forms  the  basis  for  explaining  the  spontaneous  heat  spikes  that  occur  at  low 
temperatures  when  the  atom  density  exceeds  some  critical  value. 

Thus  the  atom  concentration  obtainable  depends  on  the  initial  vacancy 
concentration.  Comparison  of  a-H2  and  C-H2  results  could  test  this  prediction.  Only  when 
B(x)  has  no  turning  point  can  the  value  of  x  exceed  the  minimum  in  the  curve  B(x)  that 
occurs  at  X  =  4/3  and  produce  an  atom  buildup  curve  like  that  seen  in  1963  by  Sharnoff 
and  Pound  [22]  and  more  recently  by  Collins  [23].  Sharnoff  and  Pound  observed  D2 
containing  1%  T2  while  Collins  studied  H2,  HD,  and  D2  containing  2%  T2.  Note  that  at  the 
critical  point,  the  saturation  value  of  x  if  there  was  no  correction  for  vacancy  filling  would 
be  0.770  rather  than  4/3.  Thus  the  saturation  level  is  higher  by  a  factor  of  1.732  =  V (3). 

The  problem  of  atom  production  in  a  solid  hydrogen  host  containing  tritium  was 
reformulated  in  terms  of  a  local,  inhomogeneous,  production  of  atoms  that  can  then 
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diffuse  into  the  bulk  of  the  solid,  forming  a  uniform  background. 

The  basic  unit  in  the  model  is  a  bubble  that  contains  all  the  atoms  created  by  a 
single  beta  decay.  After  creation  of  the  bubble,  two  effects  are  considered:  the  boundary 
expands  due  to  atom  diffusion  and  atoms  can  recombine.  There  is  no  flux  of  atoms 
across  the  bubble  boundary.  To  estimate  the  initial  size  of  the  bubble,  an  extrapolation 
of  the  measured  range  of  electrons  in  solid  hydrogen  and  an  estimate  of  the  distance 
traveled  by  20  eV  secondary  electrons  is  used.  The  size  estimate  is  that  800  atoms  are 
produced  in  a  region  occupied  by  5  x  10®  molecules. 

The  bubble,  small  in  size  initially,  contains  atom  densities  higher  than  ever 
observed  before  in  a  hydrogen  host.  The  time  constant  that  characterizes  the  lifetime  of 
the  bubbles  is  the  diffusion  time,  Tj.  In  that  time,  the  atom  density  decreases  by  a  factor 
of  approximately  10,  the  reduction  coming  from  recombination  and  expansion  caused  by 
diffusion,  with  expansion  making  the  dominant  contribution.  The  initial  high  atom 
density  results  in  spins  that  are  unobservable  by  EPR  and  negates  the  use  of  EPR 
experiments  to  determine  the  number  of  atoms  in  a  hydrogen  host.  The  spins  that  are 
observed  are  free  spins,  not  necessarily  characteristic  of  the  other  spins  in  the  solid  and 
cannot  be  used  to  determine  directly  the  population  of  spins. 

Through  analysis  of  the  bubble  geometry,  a  phenomenological  theory  that  explains 
and  brings  together  various  pieces  of  data  has  been  provided.  This  includes  the  materials 
dependence  of  the  data  of  Collins,  et  al.  [14]  and  addresses  the  concept  of  the  missing 
spins.  It  also  gives  a  framework  and  mathematical  model  that  one  can  push  further. 

There  are  many  complications  that  were  ignored  that  are  only  appropriate  for  a  much 
more  sophisticated  model.  This  includes  the  exact  nature  of  the  damage  done  by  the 
atom  formation,  the  heat,  and  possible  nonequilibrium  density  of  defects.  In  addition,  the 
formulations  of  the  missing  spins  arguments  always  start  from  initial  random 
distributions  of  spins.  This  may  not  be  justified  at  all. 

Atom  densities  on  the  order  of  2%  seem  possible  if  the  atoms  can  be  produced 
uniformly  over  the  sample  volume.  A  significant  reduction  in  the  recombination 
coefficient  resulting  from  increased  pressure,  vacancy  filling  by  atoms,  isotopically 
substituted  and  disordered  samples,  or  a  large  applied  magnetic  field  could  result  in  even 
higher  steady-state  atom  densities. 

In  work  on  samples  confined  to  two  dimensions,  estimates  of  the  molecular 
diffusion  coefficient  can  be  made.  The  data  is  in  the  fast  diffusion  limit  and  the  hopping 
frequencies  can  be  determined  from  the  NMR  line  narrowing  but  not  the  ortho-para 
conversion  rates.  The  data  is  consistent  with  the  HD  data  and  each  set  yields  diffusion 
coefficients  of  order  D  cmVs. 

There  is  no  explanation  for  the  increase  in  the  conversion  rate  observed  as  the 
samples  age  at  fixed  temperature.  Other  workers  attributed  the  increase  they  observed  to 
clustering  of  the  ortho  molecules  driven  by  the  electric  quadrupole-quadrupole  interaction. 
While  this  might  be  true  for  their  experiments,  it  is  not  a  likely  source  of  increased 
conversion  when  the  conversion  is  catalyzed  by  unpaired  electron  spins.  From  the  NMR 
lineshape  narrowing,  there  is  no  indication  of  increased  diffusion  as  the  samples  age; 
rather  the  opposite  is  true.  The  frequency  F  appears  to  decrease  (and  hence  D)  with 
decreasing  ortho  concentration  (Xj)  or  increasing  age. 
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